and other 


Aids Computation 


A Quarterly Journal 


Edited by 
J. H. BIGELOW A. H. TAUB 
Cc. C. CRAIG J. TODD 
A. ERDELYI C. B. TOMPKINS, Chairman 
IX 
nos, 49-52 
1955 
Published by 


THE NATIONAL RESEARCH COUNCIL 
Washington, D. C. 


\ 


Math. -Econ. 
Ubrary 


@A 
47 


M 
V.9 


N 


January 1955 univenMawnber 49 


OF MICHIGAN 


JUN 14 1955 


MATH. ECON. 
LIBRARY. 


and other 


Aids Computation 


lume LX 


Published Quarterly 
by the 
National Academy of Sciences—National Research Council 


ot 
ib. 
* 
“OUNDED 18 


Published quarterly in January, April, July, and October by the National Academy of Sciences— 
National Research Council, Prince and Lemon Sts., Lancaster, Pa., and Washington, D. C. 


Entered as second-class matter July 29, 1943, at the post office at Lancaster, Pennsylvania, under 
the Act of August 24, 1912. 


Editorial Committee 
Division of Mathematics 
National Academy of Sciences—National Research Council 
Washington, D. C. 
C. B. Tompxins, Chairman, University of California, Los Angeles, California 
J. H. Biaetow, The Institute for Advanced Study, Princeton, New Jersey 
C. C. Craic, University of Michigan, Ann Arbor, Michigan 
A. Erpéuy1, California Institute of Technology, Pasadena, California 
A. H. Tavs, University of Illinois, Urbana, Illinois 
J. Topp, National Bureau of Standards, Washington, D.C. 


Information to Subscribers 


The journal is now published quarterly in one volume per year with issues 
numbered serially since Volume I, Number 1. Subscriptions are $5.00 per year, 
single copies $1.50. Back issues are available as follows: 

Volume I (Nos. 1-12, 1943-1945) $12.00 (sold only as a complete volume) 


Volume I (No. 7, 1944) ‘‘A Guide to Tables of Bessel Functions,” by H. ¥@ : 
Bateman and R. C. Archibald, 104 pp., $2.00 


Volume II and III (Nos. 13-28, 1946-1949) $4.00 per year, $1.25 per issue 
Volume IV and following, $5.00 per year, $1.50 per issue 


All payments are to be made to the National Academy of Sciences and for- 
warded to the Publications Office, 2101 Constitution Avenue, Washington, D. C. 


Agents for Great Britain and Ireland: Scientific Computing Service, Ltd., 
23 Bedford Square, London W.C.1 


Information to Contributors 


All contributions intended for publication in Mathematical Tables and Other 
Aids to Computation and all books for review should be addressed to C. B. 
Tompkins, Department of Mathematics, University of California, Los Angeles 24, 
California. The author should mention the name of an appropriate editor for 
his paper if this is convenient. 


= 


i’ 


for sec 
lus 
2. | 

matter 


itie 
itie 


i 
i 


1 

abilitie 
inc 


ab 
ab 


q 


Coupon Collector’s Test for Random Digits 


1. Introduction. Increasing use of random numbers, especially in Monte Carlo 
procedures and in large computing installations, has served to focus attention on 
the various tests for randomness. KENDALL and BABINGTON-SmiITH#! list four tests 
for so-called local randomness. While not giving the coupon collector’s test (to 
be described below) a place in their now classical list of four tests, they did use 
a modified coupon collector’s test in same of their investigations. 

In an ordered set of digits, say, one may count the length of a sequence, 
beginning at a specified position, necessary to give or include the complete set 
of all ten digits. Or one may count the length required to give a set of k, k < 10, 
different digits. The distribution of these observed lengths for different initial 
positions can then be compared with a theoretically computed distribution. Such 
a test will be called the coupon collector’s test from an analogy with certain sales 
promotion schemes. 

The theoretical distribution may be computed from formulas given by H. von 
SCHELLING,? which formulas hold for the case where the individual category prob- 
abilities might be unequal. For the random digital case with category prob- 
abilities equal to 1/10, von SCHELLING’s formulas simplify readily and may be 
conveniently related to the ‘‘differences of zero.’’ These latter quantities are tabu- 
lated by FISHER and YATEs? up to sequences of length 26. But in using the coupon 
collector’s test for a complete set of all 10 digits it has been found that the mean 
of the length distribution is slightly greater than 29, and a table of probabilities 
associated with the sequence lengths 10 to 26 inclusive would hardly give a realistic 
picture of the entire distribution. 

The present author has therefore extended this tabulation, and exact prob- 
abilities are given for sequence lengths 10 < m < 35 and approximate probabilities 
for sequence lengths 36 < m < 75. The probabilities were computed from the 
relation 


@- 


~ & 
and are listed in Table 3. 

If one were interested in sequence lengths necessary to obtain five different 
digits, the mean of this distribution is approximately 6.46. The range 5 to 26 
inclusive available from FISHER and YaTEs* might be sufficient here. 

2. Sequence Lengths for Decimal Expansions of x and e. It would be a simple 
matter to program a large digital computing machine so that it would tabulate 
the distribution of the sequence’ lengths needed for complete sets for a given 
ordered digital collection. However, the author did not have such a digital com- 
puting machine available, and he made a tabulation by hand for the decimal 
expansion of x. The 2035 decimal approximation to x given by GEORGE W. 
REITWIESNER‘ was used as the raw material for this count. Beginning with the 
initial position 3 in x ~ 3.14159---, it was recorded that a sequence length of 33 
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positions was needed to get all the ten digits. Beginning anew with the thirty- 
fourth position digit (which is a 2), it was recorded that a sequence of 18 positions 
was needed to get a complete set of all the 10 digits. Continuing this procedure, 
67 sequences of complete sets were obtained, plus an incomplete sequence (at 
the end of the decimal expansion) of length 15. It was considered advisable to 
make the sequences non-overlapping as described above since there is considerable 
dependence among the set of sequence lengths if every position in the decimal 
expansion of is regarded as a new starting point. 

The sequence lengths for x are also included in Table 3. This tabulation 
for + was checked by Mr. Wayne Jones of the Department of Defense, Wash- 
ington, D. C. 

Mr. JoNEs also made a tabulation based on the decimal expansion of e. 
REITWIESNER‘ gave a 2010 decimal approximation to e. An additional 490 places 
was given by METROPOLIS, REITWIESNER and von NEUMANN.’ Mr. Jones found 
82 complete sequences using 2486 digits in the expansion of e. This tabulation is 
also given in Table 3. The author desires to thank Mr. Jones for this count. 

3. Statistical Tests. The mean and the standard deviation of the theoretical 
distribution may be computed from results given by von SCHELLING? or FELLER.® 
These theoretical values and the corresponding observed values for x and e are 
given below. 


TABLE 1 
Observed 
Theoretical 
Mean 29.29 30.16 30.32 
Standard deviation 11.21 11.83 10.64 


To use a chi-square test, it is desirable that the expected values all exceed 10 
in size. Since the sample size for is small (67) some grouping of the sequence 
lengths is necessary to meet this desired minimum. The following results were 
obtained for a convenient grouping. 


TABLE 2 
e 
Sequence lengths, n Observed Expected Observed Expected 

10-19 13 11.604 12 14.202 
20-23 13 11.720 11 14.344 
24-27 9 11.491 14 14.064 
28-32 5 11.480 15 14.050 
33-39 13 10.195 17 12.477 
40 and over 14 10.510 13 12.863 
Totals 67 67.000 82 82.000 
Chi-squared test 6.436 2.826 


values 


Neither of these chi-square test values is unusually out of line. It has been 
previously reported®:’ that (using a sample of 2000 digits for e) excessive flatness 
in the single frequencies was noted, and an indication was obtained that the single 
digits in e are “‘non-random.”’* Apparently, this phenomenon did not reflect itself 
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ty- TABLE 3 
nn Table of Probabiliti and Empiri 1 Distrit ti 
(at Observed for 
to n Pn 
ble 10 .0003 6288 0 0 
nal 11 .0016 3296 0 0 
12 .0041 9126 4 0 0 
13 .0080 9315 2 0 0 
14 .0130 4560 8576 0 1 
15 .0186 3435 9744 1 2 
| 16 .0243 5958 6451 2 1 1 
ces 17 .0297 8461 8864 6 2 
ind 18  .0345 7819 0373 1264 3 2 
wie 19  .0385 2892 7611 5744 2 4 
bt! 20 .0415 3577 5577 4998 4 2 3 
- 21 0435 8654 2461 1780 8 2 2 
22  .0447 3311 6259 6932 2752 3 3 
are 23  .0450 6836 4358 6388 8896 6 3 
24  .0447 0706 5704 2485 9072 2 3 
25 .0437 7151 9771 4451 888 1 1 
26 .0423 8153 3617 2618 0440 544 3 5 
27 + 0406 4806 4094 7299 5986 8 3 5 
28  .0386 6968 0608 0430 0677 8256 2 3 
29  .0365 3106 7596 0890 3842 8456 2 1 
fi 30 .0343 0291 1584 0099 4076 1298 5728 0 4 
ace 31 .0320 4266 1164 2497 4751 5573 3056 1 2 
ere 32  .0297 9578 2029 8315 1051 7926 5414 4 0 5 
33  .0275 9724 1577 4565 5030 9198 2083 2 1 2 
34 .0254 7304 5949 3494 9321 3424 0393 4016 2 2 
35  .0234 4171 8456 6112 5667 0619 1553 5264 2 3 
36  .0215 1565 5696 6141 0012 3 3 
37 0197 0233 0293 7275 2140 2 5 
38 0180 0533 0690 9430 5978 3 1 
39  .0164 2524 1844 5333 7918 0 1 
40  .0149 6037 8429 7183 4300 1 3 
41 .0136 0738 6073 5944 1433 1 1 
42  .0123 6172 7525 9630 8189 1 1 
43 .0112 1807 0507 6953 1223 6 1 
44 0101 7059 2895 9431 7444 1 0 
45  .0092 1321 9356 5092 1003 1 1 
46  .0083 3980 1802 1014 6739 0 0 
een 47 0075 4425 4318 8464 5255 0 1 
ness 48 .0068 2065 1566 0407 8968 0 1 
ngle 49 0061 6329 8170 0629 7386 0 0 
self 50  .0055 6677 5325 1020 3197 0 0 
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TABLE 3—Continued 


Observed for 
Pa e 
51 .0050 2596 9683 5475 6580 0 1 
52 .0045 3608 8658 5862 1077 0 0 
53 .0040 9266 5455 6666 7056 0 1 
54 .0036 9155 6480 2793 7180 0 0 
55 .0033 2893 3218 7175 0148 0 0 
56 .0030 0127 0238 7102 7949 0 0 
57 .0027 0533 0592 0946 3793 0 0 
58 .0024 3814 9607 8519 0648 0 0 
59 .0021 9701 7828 5704 8789 1 0 
60 .0019 7946 3656 1777 4941 1 0 
61 .0017 8323 6124 7283 4517 0 0 
62 .0016 0628 8101 7164 6560 0 0 
63 .0014 4676 0128 6430 1251 0 0 
64 0013 0296 5041 3524 9640 0 0 
65 .0011 7337 3456 8177 1422 0 0 
66 .0010 5660 0172 2241 9129 0 1 
67 .0009 5139 1491 6632 0338 0 0 
68 .0008 5661 3473 2828 0290 0 0 
69 .0007 7124 1073 6049 1625 0 0 
70 .0006 9434 8154 4810 4916 0 0 
71 -0006 2509 8310 7043 5014 0 0 
72 .0005 6273 6471 7289 2795 0 1 
73 0005 0658 1228 5628 1531 : 0 0 
74 .0004 5601 7836 1436 4579 0 0 
75 .0004 1049 1841 9142 4169 0 0 
76 .0036 9745 8744 5702 0432 1 (77) 0 
and over 
Total 67 82 


in materially changing the characteristics of the sequence length distribution for 
the coupon collector’s test. Some question arises as to whether the single frequency 
test and the coupon collector’s test are independent, and ale which test has the 
greater power. 

The chi-square test values in Table 2 were calculated by assuming that the 
sequence lengths for complete sets of digits are independent draws from a known 
(infinite) multinomial probability distribution. (Null hypothesis.) The alterna- 
tives would include unspecified sorts of dependency and other underlving prob- 
abilities different from those given in Table 3. 

RoBErT E. GREENWOOD 
The University of Texas 
Austin, Texas 

1M. G. Kenpatt & B. BasincTon Sirs, “Randomness and random sampling numbers,” 
Royal Stat. Soc., Jn., v. 101, 1938, p. 147-166. 

?H. von SCHELLING, “Auf der Spur des Zufalls,”” Deutsches Statisches Zentralblatt, v. 26, 1934, 


p. Sag ty Also, “Coupon collecting for unequal probabilities, ” Amer. Math. Mon., v. 61, 1954, 
p. 306-311 
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*R. A. FisHer & F. Yates, Statistical Tables for Biological Agriculture and Medical Research. 
3rd edition, London, 1948, Table XXII. 

4 GrorGE W. REITWIESNER, “An ENIAC determination of x and ¢ to more than 2000 decimal 
places,” MTAC, v. 4, 1950, p. 11-15. 

5.N. C. Metropo.is, G. REITWIESNER, & J. von NEUMANN, “Statistical treatment of values 
of first 2000 decimal digits of e and x calculated on the ENIAC,” MTAC, v. 4, 1950, p. 109-111. 

6 W. Fetter, Probability Theory and Its Applications. Volume 1, New York, 1950, p. 175-181. 

7F. GRUENBERGER, “Further statistics on the digits of ¢," MTAC, v. 6, 1952, p. 123-134. 


A Method for the Evaluation of a System 
of Boolean Algebraic Equations 


With the advent of large scale electronic devices whose logical design is de- 
scribed by a system of Boolean algebraic equations, a method to mechanize the 
evaluation of such a system and shorten this evaluation with respect to time will 
be increasingly useful. Such a method will be described in this paper. 

The problem may be described as follows: Given a set of m variables, Q*, 
(k = 1,2, ---,) each of which may take on the value 1 (true) or 0 (false) at 
any time ¢; then the value of any Q* at time ¢ + 1 may be defined by the system 
of Boolean equations 


(1) R# = 
(2) Sé = ge(Q:*) 
(3) = SA) 


' where 1 < q < n. For example, the recirculation loop of a dynamic flip-flop may 


be defined simply by 
= (04, S#) = Re. 


In another system, a more complex definition 
= R#, S#) = Of-RE-S# + + OF-Ré-S# 


may be taken, where R;* and S/ are the two inputs to flip-flop Q*. 
We shall use the symbols for conjunction, disjunction, and negation 


Q'-e “Q! and Q?” conjunction 
e+e “Q' or @” disjunction 


Q “Not Q?” negation 
which are defined by the truth tables? , 


0 0 0 0 1 

0 1 0 1 i 

1 0 0 1 0 

1 1 1 1 0. 
A “term” is defined as one or more variables conjoined together, e.g., Q'-Q*-@; 
and an “equation” as M terms, T,, (m = 1,2, ---, M) disjoined together, e.g., 


Q'-Q?-@® + Q'-Q*. Now note that the value of a term is zero if any variable in 
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that term has the value zero, and the value of an equation is one if any term in 
that equation has the value one. We will deal only with the case where all of the 
expressions on the right are in normal disjunctive form, but this represents no 
restriction on the method since any equation in the above form may be so written. 

Our method for solving this system of equations involves the use of punched 
cards. While this description is in terms of IBM punched cards, the method is, 
of course, not restricted to this type of card. Cards and a key punch are the only 
tools required. In fact, the key punch is required only to initially punch the cards 
as described below, and is not used thereafter in the actual process of evaluation. 

Two types of cards must be defined and punched in order to implement the 
evaluation. Cards of type I each represent the m term of the equation for R* 
or S* and are labeled “‘R*T,,”’ or “‘S*T,,”; these are the ‘“‘equation’’ cards. Cards 
of type II represent the value of Q* (1 or 0) and are labeled Q*T or Q*F, respec- 
tively ; these are the ‘“‘value’”’ cards and one or the other is selected for use at time 
t for each k. The problem may now be restated as: given the set of cards of type II 
containing the appropriate card Q*T or Q*F for all Q* at time ¢, and further given 
the equations as represented by the deck of cards of type I, construct a set of 
cards of type II representing the value of each Q* at time ¢ + 1. 

The twelve rows and 80 columns of the IBM card may be thought of as a 
single pair of rows and 480 ‘“‘columns.”’ Let this column number correspond to k 
and let us represent variable Q* by one of three possible combinations of punches 
or no punches in the pair of rows of column k. 

Thus cards of type I are constructed such that the card R*T,, has column g 
punched as follows: 


i) a punch in row 1 if Q* is a variable in term m of the equation for R*, 
ii) a punch in row 2 if Q* is a variable in term m of the equation for R*, or, 
iii) no punch in column q if neither Q* nor Q¢ appears as a variable in term m 
of the equation for R*. 


The card S*T,, is constructed in an entirely similar way. 

Cards of type II are constructed in a manner which is somewhat analogous 
to the above. The same 480 columns of a single pair of rows is considered. Now, 
however, there -is one and only one punch in column &, and all rows of all other 
columns are punched ; i.e., 


i) the punch in column 2 is in row 1 if Q* is true . . . in this case the card is 
labeled Q*T; we 
ii) the punch in column & is in row 2 if Q* is false (Q* is true) . . . in this 


case the card is labeled Q*F; 
iii) both rows of all other columns are punched. 


A “‘deck”’ of cards of type II consists of all of the cards Q*T and Q*F. A “‘set’’ of 
these cards at time ¢ consists of one and only one card Q*T or Q*F representing 
the value of each Q* at time ¢. Thus there are twice as many cards in a deck of 
type II as in a set of type II. It will be convenient to have two decks of type II 
so that one set of cards may be made up for time ¢ from one deck and, using this 
set and the entire deck of cards of type I, a set of cards of type II for time ¢ + 1 
may be compiled from the other deck. 
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We shall call the set of cards of type II, for time #, the ‘‘this time table” and 
those for time ¢ + 1 the “next time table.” Note that after we have formed the 
complete next time table, we may redefine it as the new this time table and proceed 
to form a new next time table for time ¢ + 2 from the (now) unused deck of type 
II. This process may continue as long as required, and is essentially that of 
evaluation of the equations. 

The process is carried out as follows: The this time table is stacked and re- 
mains undisturbed (but not unused) until the next time table is completed. The 
entire this time table is physically aligned so that each position on the card 
representing each row and column is directly in line with the corresponding 
position of all other cards in the set. With these cards so arranged, it is easy to 
determine which positions are punched in every card of the set ; this is most simply 
done visually by simply sighting through the entire deck. Further, we may now 
observe that we can sight through one and only one of the two rows in every 
column. That this is the this time table for time ¢ has the significance that if we 
can sight thru row 1 of column k, then Q/ = 1, or, if we can sight thru row 2 
of column k, then Q¢ = 0. Thus by simply aligning the cards of the this time 
table, we can read the values of Q/ required to evaluate ¢ (Eq. 3). 

We may now evaluate R/# (Eq. 1) by the following measures if the cards 
R'T,, of type I are arranged in ascending order of k and m: Place card R'T, in 
alignment with the this time table as above. If we can sight thru the deck of 
cards thus formed wherever R'T; is punched, then R? = 1 and it is not necessary 
to examine cards R'T>, R'T;, ---, R'Ty, at this time. If not, we remove R'7; 
and proceed to examine R'T; in the same way. This process continues until either 
1) we find R'T,, satisfying the condition that we can sight thru the deck con- 
sisting of it and the this time table, in which case R? = 1, or 2) we have examined 
all R'T,, without satisfying the condition, in which case R? = 0. Having thus 
evaluated R/, we proceed in the same manner to evaluate S;. When this is done, 
we are now in a position to solve for Q',,: from equation 3. In exactly the same way 
we proceed to evaluate in order R?, S?, S?, Pur, ---, Re, Se”, QP 
The values of the Q’s are recorded as the next time table by forming the proper 
set from the available deck of cards of type II. The values of R and S need not 
be retained. 

The above description in which IBM cards are employed requires that 
N < 480. It is easy to see that this number may be increased by increasing the 
number of cards of type I and type II appropriately. However, the more usual 
problem is to reduce the number of cards employed when 1 is sufficiently less 
than 480, this being the more usual case. 

The number of cards of type II employed is perfectly straightforward: If 
n < 480, then the number in a deck is 2m and that in a set is m; if 480 < m < 960, 
then the number in a deck is 4m and that in a set is 2”; and so forth. But if is 
sufficiently small, e.g., m = 240, then we may make certain types of combinations 
of the equation cards (type I). For example, we may represent R*T,, and S*T,, on 
the same card, or we may represent two terms of R* on one card. Such combina- 
tions of cards serve to reduce the time required to carry out the evaluation, 
where 7 is sufficiently small, without loss of generalization. They require certain 
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obvious modifications of the configuration of the value cards (type II) in the 
form of duplication of information. 

When the problem actually is solved as described above, it is convenient but 
not necessary to have such devices as a box of some transparent material in which 
to stack the this time table, perhaps with a strategically placed light source to 
assist in observations. In fact, the design has been completed for'a device to 
completely mechanize the process described so that the human operator has only 
to feed cards to a card reader with the solution carried out and the results re- 
corded automatically. 


J. A. PostLey 
Hughes Research and Development Laboratories 
Culver City, California 


1 Hans REICHENBACH, Elements of Symbolic Logic. The Macmillan Company, New York, 1947. 
bi as used in this paper is that of Reichenbach except that his symbol ‘‘v” is replaced by 

“ +”, 

2 E. C. NEtson, “An algebraic theory for use in digital computer design,” Transactions of the 
IRE Professional Group on Blectronic Computers, vol. EC-3, p. 12-21, September 1954. 

3 ENGINEERING RESEARCH AssociATES, High-Speed Computing Devices. McGraw-Hill Book 
Company, New York, 1950. 

* ArTHuR W. Burks, Don W. WarREN, & JESSE B. Wricurt, “An analysis of a logical machine 
using parenthesis-free notation,’”’ M7‘AC, v. 8, 1954, p. 53-57. 


Tables for the Determination of Fundamental 
Solutions of Equations in the Theory of 
Compressible Fluids 


1. The pseudo-logarithmic plane. It has been shown'** that it is useful, 
when applying the hodograph method,‘ to consider the stream function y and 
the potential function ¢ of compressible fluids in the so-called pseudo-logarithmic 
plane. In the case when the pressure density equation is = ap’, o, y being con- 
stants, and for subsonic flows, the cartesian coordinates of the pseudo-logarithmic 
plane are 


and @ where M is the Mach number and @ is the angle which the velocity vector 
forms with the positive x direction of the physical plane (t.e., the plane in which 
the flow actually takes place).® 


The equations for the potential (¢) and the stream function (W) assume, when 
considering the flow in the X, @-plane, the form 


(2) grr + + Jo = 0, 


(3) da + vo + (F(A) = 0. 


(12) 
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Here 1 = 1(A) = p~ (1 — M”) is a known function* of \. If we introduce so-called 
modified stream and potential functions 


(4) 

(5) = Kd, 

(6) c= (1 — + — 
then we obtain for ¥* and ¢* equations 

(7) + vo + 4Fi(A)y* = 0, 

(8) + + = 0 

where’ 


(9) F, = (—1) 


64 (1 — M?)* 
= (37-1), a=7-3. 


2. Singular solutions of type S. By a fundamental solution of a linear differ- 
ential equation we mean a function $*(A, 6; Ao, 80) which for a fixed (Ao, 40) is a 
solution of the differential equation and which possesses a logarithmic singularity 
in the vicinity of (Ao, 40). In the case of equations (2) and (3) singularities of th's 
type (denoted as singularities of type $)* can be obtained by continuing the 
coefficients F,; and F, to complex values of \ and then forming for ¥* the expression 
(10) $A[log (A — ro)? + (6 — + B 


where 


10(H"A) 10(H-"A) 


af t f=Z-Z, 


Z=rA+ 0, Z=r—, Zo + io, = ro — Ho. 


We obtain the corresponding singularity ¢* replacing in A and B functions H 
and F, by and F,, respectively.® 

In order to evaluate numerically the integrals in (11) and (12) and correspond- 
ing expressions for ¢* we need the tables of Re(F,), Im(F,), x = 1, 2 for complex 
values of the arguments \. 
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ImQ) Re(T)  Im(7) Re(F;) Im(F;)  Re(F:) Re(Q 
a) 0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 —0.34 
—1.1513 0.0000 .9750 0.0000 .0017 0.0000 — .0068 .0000 


RSE 

233 
338 
orn 
£38 

2 

3338 


2617 9788 — .0130 0015 — .0026 — .0061 
3271 -9808 — .0157 0003 0016 — .0008 — .0064 
3925 -9833 — .0180 0003 0015 + .0010 — .0062 
4579 -9859 — .0200 — .0007 0014 0026 — .0056 
5233 -9888 — .0216 — .0011 0011 0040 — .0045 
5887 9918 — .0228 — .0012 .0008 0049 — .0032 
6541 -9949 — .0235 — .0014 0005 0054 — .0017 
7195 -9980 — .0238 — .0014 0001 0055 — .0003 
7849 -9989 — .0237 — 0013 — .0001 0054 0001 
—0.8047 -0000 -9473 .0087 — .0348 
0654 -9478 — .0076 0081 0027 — .0329 — .0112 
-1309 -9496 — 0151 0068 0052 — .0270 — .0210 —0.21 
1963 -9526 — .0219 0046 0069 — .0185 — .0276 =" 2 
2617 -9565 — .0282 0022 0077 — .0087 — .0308 - 2 
3271 -9613 — .0339 — .0003 0077 0013 — .0307 =— 2 
3925 -9667 — .0385 — .0024 0069 0098 — .0276 
4579 -9726 — .0422 — .0041 0056 0164 — .0223 vocal 
5233 — .0450 — .0053 0208 — 0158 
5887 9852 — .0469 — .0059 0023 0231 — .0091 
6541 9916 — .0479 — .0058 0006 0233 — .0026 
7195 9981 — .0479 — 0055 — .0009 0218 0034 
7849 1.0042 — .0473 — .0049 — .0021 0193 0083 
—0.6020 -9161 -0000 -0000 — -0000 
0654 9172 — .0129 0238 0095 — .0961 — .0383 
-1309 9206 —, .0252 0181 0169 — .0732 — .0684 
-1963 -9261 — .0366 0212 — .0416 — .0854 
-2617 -9332 — .0466 0025 0221 — .0092 — .0885 — 
3271 -9416 — 0551 — .0045 -0200 0190 — .0802 “4 
3925 -9510 — .0618 — .0098 -0162 0396 — .0643 > 
4579 -9608 — .0669 — .0130 0115 0523 — .0454 ¥ 
$233 — .0703 — .0143 0575 — .0262 
5887 -9810 — .0721 — .0143 0022 0570 — .0088 - 4 
6541 9909 — .0727 — .0134 — .0015 0527 0059 - 
7195 1.0004 — .0719 — 0115 — .0043 0456 0173 - 
7849 1.0094 — .0701 — .0093 0064 0371 0257 - 
—0.5281 -9000 0401 -0000 — .1621 0.0000 
—0.4582 -0000 -8810 -0000 -0629 -0000 — .2556 -0000 
-04 88181 — .0120 .06018 01654 — .2441 — .0679 
0654 -8830 — .0195 0558 0261 — .2266 — .1072 
1309 -8888 — .0380 0381 — .1521 -1818 
1963 8980 — .0545 0158 0515 — .0608 — .2081 
20 89856 — .05522 01459 05137 — .0564 — .2079 
—0.4582 -2617 -9094 — .0684 — .0041 0483 0194 — .1939 
3271 9227 — .0794 — .0184 0388 0757 — .1541 = 
3925 9368 — .0877 — .0265 0272 -1069 — .1064 
4579 9511 — .0932 — .0293 0156 -1172 — .0612 
5233 9653 — .0967 — .0289 -0060 -1146 — .0226 a 
5887 9791 — .0978 — .0260 — .0018 -1030 0072 “* 
6541 -9923 — — — .0073 -0874 0292 
7195 1.0048 — .0953 — 0176 — .0112 0701 
7849 1.0163 — .0919 — 0133 — .0135 .0529 0535 


8 
1 


Re(a) 
—0.3466 


-1780 


.1779 


.1778 
1777 
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Im(A) 
-0000 


Im(T) 


01751 


Re(F;) 


Im(Fi) 


= 


il 
0000 .1425 .0000 .0000 
04 84168 — | .12942 04434 — .1849 
0019 0654 .8436 — .0285 1195 0689 — .2877 
0038 10 84825 — .0426 .09286 .09366 — .3879 
0051 1309 4.8535 — .0547 0674 .1070 — A416 
0061 17 86179 — .06835 .03083 .11160 — 4561 
‘0064 1963 .8682 — .0769 0095 — Asli 
‘0062 2617 8860 — .0945 — .0295 0877 — 3491 
‘0056 3271 — .1075 — 0496 0596 — .2341 
‘0045 3925 — .1163 — 0556 0337 2244 — .1293 
.0032 4579 — 0133 19° 
5233 .9625 .1240 — 0468 — .0016 .1845 0088 
.0003 5887 .9799 — .1238 — 0385 — .0114 (1517 0466 
.0001 6541 9961 — .1215 — 0304 — 0177 -1189 0699 
7195 1.0110 — .1175 — 0224 — .0212 0887 0829 
0000 7849 «1.0246 .1123 ~ 2227 0889 
-0210 —0.2682  .0000 .8000 .0000 .2837 .0000 —1.181 .0000 
— .3275 .8889 — .1390 — .1050 .0718 4285 — .2722 
-0308 — .2560 .5237 .9625 — .1515 — 0643 — .0179 .2507 0743 : 
7193 1.0188 — .1388 — 0248 — .0334 
— 2559 .1968  .8375  — .1049 — .0365 .2015 1749 8151 
— 2558  .5889 — .1491 — 0991 — .0272 .1907 1075 
0091 — .0000 .7921 .0000 3217 .0000 —1,342 .0000 
0653 .79830 — .0408 2413 — 9944 — .7737 
0026 2623 ©8633 — +1255 — 0923 .1312 3900 — .5160 
0034 3926 9155 — .1472 — 0966 0270 3846 — .0952 
.7930 — .0163 .30819 08095  —1.284 —0.3456 
0000 0502 .7956 — .03165 27223 14881 —1.127 — 6331 
0383 09987 .8056 .06055 .15893 23065 — .6380 — .9700 
‘0684 14 81717 — .08116 06324 24157 — .2327 —1.002 
0854 
‘0885 — .1313 .8144 — .0768 —.3171 —1.010 
6540 1.0018  — .1450 — 0359 — .0317 .1397 1247 
0802 — 2550  .93988 .1510 — 0813 — .0010 .3184 0113 
0643 — 2549 .7852 1.0343 — .1316 — 0156 — .0330 0618 .1289 
0454 — .8863 .1563 — .1371 5543 .2191 
0262 
0088 — 1815 S791 .1737 — 0577 — .0460 2211 .1807 
— 1813 .7330 .0000 —3.307 
0059 — 1806  .7434 .0590 4632 — .5113 — 1.908 —2.230 
0173 — .1792 1.0095 — .1670 — 0372 — .0479 .1864 
0257 — 1787  — .1739 — .1785 0518 174 — .1713 7 
0000 — 1785  .2624  .8434 #— .1620 .0022 8638 — .5537 
0000 — 1783 0000 .7223 .0000 .78138 —3.862 | 
0679 0242 .7339 — .02275 .72795 24278 —3.100 — 1,068 
1072 04 73691 — .03669 64325 36818  —2.718 —1.614 
1818 05006 .7390 — .0462 .57508 44169 —2.412 —1.929 
10 75633 — .08624 19955 56744 —.7558 —2.422 
10026 .7569 — .0868 .19337 56525 —.7292 —2.411 
2192 .8198 — — .20073 8695 — .9593 
1939 5893-9886. — ,1738 — 0544 — .0474 1855 
1541 7197. 1.0279 — .1587 — 0236 — .0465 1792 
- 7856 1.0449 — .1493 — 0128 — .0434 0523 
«3921S 9097S — «1800 — «1414 — .0017 5529 0283 
0072 4576 — .1813 — 1064 — .0293 4096 
0292 1968  .8072 — .1414 — .1816 .3108 ‘8132 —1.229 
0443 1307 — .1074 0116 5161 0452 —2.161 
0835 5244 9655 — .1792 — 0771 — .0426 .2951 .1722 


0236 


-1739 


Re(7) 


Im(T) 


-1389 


-2100 
-1875 
-2031 


Re(F;) 


— .1904 
1.2236 
2.1870 

— .2491 

— .0347 

— .3484 


& 


oo 


— .2005 
—2.3558 


Im(Fi) 
4119 


a 
a 


21.9170 


8.07236 
4.41078 
92169 
— .1587 
— .0276 
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Im(\) 


Re(F:) Im(F:) 


8786 — 1.659 
—5.309 
—9.677 0000 

9739 0952 

1333 -2470 
1.416 — .2377 
—2.802 —7.041 
—6.153 —6.487 
—5.371 —6.113 
—3.538 —7.002 
— .7528 —6.596 
2.214 — 2.382 
6739 2390 
2025 -2650 

0785 -2281 

3015 -2850 

4645 -2938 

2.058 — 1.160 
2.099 —3.769 
0344 -2065 
—6.479 
—21.26 d 
—3.409 —23.27 
1.098 4646 
-1134 
— 39.90 -0000 
-1732 3337 
1.694 4255 
5.036 —22.47 
-7184 -5060 
— 16.07 — 36.92 
—4.884 —33.74 
2.574 — 28.37 
0567 
2627 
3.215 -2607 
7.586 —2.639 
-4342 4732 
0085 2414 
—88.02 -0000 
1.111 -6361 
— 103.8 -0000 
1028 -3228 
1513 -3584 
—26.78 —88.88 
1.632 -7649 
28.27 —32.98 
9.923 — .8918 
—255.0 .0000 
—250.3 —39.89 
—101.9 — 169.4 
26.47 —105.1 
36.54 —35.04 
28.88 —17.62 
13.68 —2.415 
-7174 6559 
10.15 1.591 


Re 


—0. 


12 
— 16444 .7878 
— 1473 .0000 .7000 .0000 .0000 t 
— 1160 .6570 .0000 .0000 
— - — .0087 
— 1136 1.0176 - — .0645 
— 1127 .2642 .8304 - 
— 1123 0629 .6762 — .0867 .7043 1.5800 
— .1116 04073 .66325 — .0589 1.43653 1.4304 re 
05 66858 — .06762 1.17971 1.48423 iZ 
06 67341 — .08015 86799 1.56464 
0798 — .1062 25015 1.50623 
15819  .7496 — .1676 — 60628 
J 
— .1112 .3902 .9072 — .2134 
— 1108  .5883  .9964 — .1963 
1.0376 — .1779 
— 1101 — .2056 
— 1094 .1982 .7807 — .1874 
— .7243 — .1516 9191 
ee — 1087 .7854 1.0564 — .1664 — .0536 
— 1061 .1216 1.4992 = 
— .0807 .0000  .6000 .0000 4668 .0000 
— .0628 .0468  .5909 — .1042 1.0202 5.0430 
— 0613 .3248  .8680 — .2443 — .2913 — .0969 
— 0610 .6497 1.0255 — .2056 — — .0791 
— .0000 .5549 .0000 8.5619 .0000 
— .0590 .5957 1.0062 — .2159 — — .0879 0. 
— 0572 .2700 — .2454 — 4358 — .0673 
— 0554 .5913 — — .7605 4.9886 
— .0528 .3838  .9062 — .2472 — — .1203 
— 0527 .03085 .5521 — .0828 3.78079 7.85575 
04 — .1030 1.4420 7.2060 
04936 .5756 = .1213 — .16725 6.18058 
as — 0506 .7115 1.0471 — .1967 — 0138 — .0723 
— 0499 .9821 — .2311 — 0719 — .1034 
— 0460 .2024 .7650 — .2459 — .7883 
— 0457 .1260 .6769 — .2170 — 1.6765 .8002 
— 0444 4555 .2442 — 1203 — .1201 
Bh — 0439 .7850 1.0692 = .1831 — .0006 — .0629 
— .0423 0000  .5000 .0000 .1836 .0000 
oc — 0407 .3229 .8670 — .2586 — 3007 — .1398 
— .0000  .4889 .0000 21.5359 .0000 
— .0385  .6966 1.0287 — .2137 — — .0854 
— 0365 .6001 1.0122 — .2236 — — .0950 
— 0349 .4924 — .0868 .1402 .1818 
— 0338 .2757 .8307 — .2648 — 4330 — .1516 
— 0301  .0487 — .1655 —5.3652 7.5268 
— 0275 1201 .6578 — .2394 —2.2320 .4906 
— 0256 .0000 .4309 .0000 51.3924 .0000 
00192 .4312 — .0098 50.4590 7.79669 
01358  .4438 — .0655 21.74356  33.62862 Sta 
02930  .4773 — .1227 —4,06835 Sta 
05023 .5265 — .1714 —7.00589 
06358 .5557 — .1926 —5.77547 Off 
10204  .6280 — .2309 —2.98574 
— 0254 .3774  .9047 — .2647 v. 
— 4.1224 ~=—«.6546 — .2622 
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Re(a) Im(Q) Im(T) 


— .0231 0691 — .2052 
— .0220 -7055 = 1.0514 — .2067 
— .0201 0000 -4000 
— .0186 -1943 -7518 — .2736 
— .0172 5367 9901 — .2436 
— .0168 0108 3939 — .0695 
— .0151 0367 — .1710 
— 0114 4538 9511 — .2617 
— .0103 -7851 1.0768 — .1919 
— .0072 0175 3759 — .1368 
00225 .0000 -1988 0000 
07135 .5483 — .2550 

12327 .6520 — .28075 

19557 .7516 — .29248 

26137 — .29218 

32557 ~—.8718 — .28652 

-38922 .9153 — .27770 

45699 .9544 — .26600 

—0.00225 .52564 .9881 — .25254 
58995 1.01533 .23890 

64888 1.03721 — .22580 

72008 1.060335 — .20940 

78546 1.07880 — .19395 

0401 — .22137 

0060 -2600 — .10704 

0.0000 .0000 0000 0000 
0248 -38623 — .2030 

0751 5435 — .25965 

-1238 6522 — .2846 

* .1967 7528 — .29514 

2612 8190 — .29424 

-3260 8724 — .28813 

9156 — .27913 

4622 95759 — .2662 

5251 9883 — .25368 

5899 1.0158 — .2398 

6491 1.0378 — .22653 

7193 1.0608 — .21012 

7857 1.0794 1945 

0019 0187 2000 


3 


0163 5196 

-3078 
0366 -1012 6062 3500 
0722 1447 6928 


3 
3 


Stanford University 
Stanford, California 


Re(F;) 
— .5637 


14 
oo 
x 


Im(F;) 
3.2184 


1.16194 
549.7400 


— 1.7443 
84843 


~ 
w 


Re(F:) Im(F;) 


27.62 —11.91 
0441 -2980 
425.2 -0000 
3.651 1.261 
-2228 4440 
165.4 —404.9 
85.83 —63.54 
5708 
— .0086 2580 
328.2 —339.0 
37290. 0000 
32.28 6.548 
9.866 3.587 
3.467 1.948 
1.733 
9555 
-6056 7419 
-2232 -4829 
1327 -4081 
0736 -3540 
0242 
— .0132 -2616 
3 6. 
3535. --2818. 
292.2 38.01 
33.23 8.931 
9.684 3.956 
3.389 2.023 
1.719 1. 
9775 9714 
5998 -7523 
4153 5435 
.2212 4875 
-1304 4107 
0716 3558 
0204 -3029 
— .0146 -2626 
285. 93.65 
112.0 62.84 
29.15 27.79 
1.200 1.362 
7.405 13.43 
1.066 6.842 
-1881 1.053 
— .7788 3.576 
— 1.200 1.871 


The eating Shin an per and the research which it reports were done under contract with the 
Office of Naval R 
1S. Bercman, “A formula for the stream function of certain flows,” Nat. Acad. Sci., Proc., 
v. 29, 1943, p. 276-281. 
2's. BERGMAN, “On two-dimensional flows of compressible fluids,” Nat. Adv. Com. for Aero- 
naulics, Technical Note No. 972, 1945. 


‘0000 
0000 
— 0103 — .0790 
844.3000 .0000 
— 9256 — .1897 | 
| — 0622 — .1173 
041 3536 7888 
487 — 11631 1449 | 
113 — 1125 — .1477 
‘002 0043 — .0674 
— .5972 .7007 
6885.0000 0000 
.2390 —7.03253 — .39862 
2650 —2.36513 — .46519 
2281 — 91111 — .35304 
2850 
2038 — .27753 
= “2472 
— 15368 | 
2065 — .12835 
479 
— :10919 
— 109458 
— 107979 
— .06838 
:2999 —21.8887 
‘0000 
.3337 
4255 —7.2871 
—2.3708 
92 — .89654 
74 — 47812 
— :27990 
37 — .17371 
2743 
3958 — 10146 — .15322 
2607 — 106264 — .12966 
639 — 03561 — 
— 101829 — .09509 
4732 — 00360 — .08023 
2414 
00605 — .06863 
— 9771 — .2243 
— 12369 — 9431 
0000 — (7018 —4.6774 
2553. — 308 
3584 —2.2466  —2.5925 
88 — 6400  —1.5496 
— 0772 — 2886 
98 — 0350 — .9719 
| 2100 .6235 
891 
STEFAN BERGMAN 
4 
04 | 
62 
A15 
6559 
591 
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*S. BERGMAN, “Two-dimensional subsonic flows of a compressible fluid and their singularities,”’ 
Amer. Math. Soc., Trans., v. 62, 1947, , P 452-498. 

*S. A. CHAPLYGIN, “On gas jets,” Scientific Memoirs, Moscow Univ., Math. Phys. Sec., 21 
(1902), p. 1-21. (English translation published by Nat. ‘Adv. Com. for ‘Aeronautics, Technical 
Note No. 1064, 1944, ard also by Brown y manda 1944.) 

5 See page 462 of the reference of note 3 

® See the reference in note 2 and equations (2.8) and (2. 5 of the reference in note 3. 

bares where symbols F and P were used for the F; and F: 
present 
8 See page G72 o of the reference of note 3. 

® We note that when we are carrying out the integrations in (11) and (12), \» are continued 
to complex values Xn + iAn, and Zn = An + iAn and Z, = dq — iA, become two independent vari- 
ables. Also see page 473 of the reference in note 3. 


Analysis of Problem Codes on the Maniac 


The Los Alamos computer, the MANIAC, has solved problems of wide variety 
during the last few years. The mathematical structure of such problems has 
ranged from differential equations (particularly partial differential equations), 
integral equations, stochastic processes, purely algebraic problems to some in 
the domain of mathematical logics. 

There are several reasons why a frequency analysis of the computer as used in 
several typical problems might be useful. From such a study one may learn of 
significant variations in such distributions from one type of problem to another. 
Further, one may reach conclusions about the selection of the computer vocabu- 
lary. Most importantly, however, one may use the quantitative results as guiding 
principles in the design of a new computer. The economy of computer design is 
connected with the question, ‘‘Is the desirability of a particular order! ,commen- 
surate with the associated electronic hardware?’’ Finally, a frequency analysis 
enables one to form accurate estimates for the “running time” of a problem; this 
information aids considerably in efficient scheduling of computer time. A knowl- 
edge of operation times for subroutines enables one to make rather good time 
estimates of lengthy problems during the forraulation stage. 

These frequency distributions can, of course, be gathered by hand. A more 
obvious way is to have the computer itself perform the analyses. A routine has 
been developed for this purpose and is called the ‘‘Code Analyzer.” 

The Code Analyzer gives the following information about a computer problem: 


1. the frequency of occurrence of each order as it appears in the code—more 
briefly, a ‘“‘static’’ count 

2. the distribution in per cent of these static counts 

3. the frequency of performance of each order during the running of a repre- 
sentative cycle of the problem—more briefly, a ‘‘dynamic’’ count 

4. the distribution in per cent of these dynamic counts 

5. the total time consumed performing each order of the vocabulary; i.e., 
(3) multiplied by the time needed to perform one such order 

6. the per cent of the total time used by each order 

7. the totals for (1), (3), (5). 


The count in (1) is obtained by simply scanning linearly through the code of 
the problem and recording the occurrence of each order. The distribution (2) is 
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the result of normalizing the distribution (1). To obtain (3), frequency of dynamic 
count, the problem must be run by means of an interpretive routine’ (in our case 
it has actually the same vocabulary as the computer) which tallies one for the 
order as each instruction proper is performed. Distribution (4) is simply a nor- 
malization of (3). The operation times for the various individual orders are stored 
as constants within the Code Analyzer, and therefore only a simple multiplication 
is needed to obtain (5) in all cases except for the shift orders which are, of course, 
of variable duration. Each shift order is examined to determine the number of 
places shifted and the time computed accurately according to a simple linear 
expression. 

Operation times were obtained by taking the time difference of a large number 
of passes through a control loop with and without each order. Dividing this time 
difference by the total number of passes, one obtains a realistic measure of the 
individual operation times that includes the access time for each instruction. 

When the above statistical information has been gathered, the percentages 
and totals are calculated and all six columns (1 to 6, above) are printed, preceded 
by a column which lists the MANIAC’s thirty-six orders. Totals (7) are printed 
below the appropriate columns. A brief explanation of the vocabulary symbols is 
given in Table 1. 

Conceptually, the idea of the Code Analyzer is very simple. In practice, this 
would also be true provided a computer had only one medium of storage which 
was sufficiently capacious to contain simultaneously the Code Analyzer and the 
problem at hand. However, this condition does not obtain for the MANIAC 
which has one thousand and twenty-four locations of electrostatic storage 
(primary) and ten thousand locations of magnetic drum storage (secondary). 

The communication between the two media gives rise to considerable compli- 
cations during the interpretative process. As a result, the code for the Code 
Analyzer itself is by no means short; indeed, it has one thousand instructions, 
apart from numerical storage. 

The procedure is as follows: First, the problem to be analyzed is loaded into 
the storage normally, as if it were going to be run. Then a key instruction is 
inserted at the end of a typical cycle to indicate the terminal point of interpreta- 
tion. With very short problems, or with subroutines, this stop is the end of the 
problem or the exit from the subroutine. The code of the problem residing in the 
electrostatic memory is then transferred to the drum at a place that the Code 
Analyzer will use as a simulated one thousand twenty-four word memory during 
interpretation. In the interest of speed, the numerical storage, both constant and 
dynamic, remains in the true electrostatic memory at all times. This is possible 
if there are several Code Analyzers, each occupying a different section of the 
electrostatic memory, so one may be chosen which does not conflict with the nu- 
merical storage of the problem being analyzed. The next step is to load the selected 
Code Analyzer into the memory. To allow for problems with prodigious storage, 
all of the code for the Code Analyzer does not reside in the electrostatic memory 
at one time. The Code Analyzer divides naturally into two parts since the count- 
ing is independent of the totaling, computing of percentages, and printing of the 
output. Thus, only the code for the counting is loaded immediately into electro- 
static storage. The second part of the code is temporarily stored on the magnetic 


1 
j 
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drum. Then, with the help of an auxiliary read tape which contains the initial 
and final addresses of the problem code, the address where interpretation is to 
begin, and the initial and final addresses of the problem storage, the Code Analyzer 
proceeds with its work and prints the results as described above. In addition 
both the dynamic and static counts for each analysis are punched on paper tape. 


TABLE 1 


A Brief Description of the MANIAC Vocabulary Symbols 


The letters are tetrads of binary digits. A, B, ---, F correspond to the hexa- 
decimal characters 10, 11, - --, 15. Associated with the pair of vocabulary symbols 
for each order are three tetrads specifying the memory location of the operand. 
There are two instructions per word; the occurrence of some orders in duplicate 
is necessary for reference to the left or right portions of the word. 


Order Vocabulary 
No. Symbol 


Add number (N) in memory (M) to cleared ac- 
cumulator (A) 

Subtract N in M from cleared A 

Recall N on magnetic tape to quotient register (Q) 

Record N in Q on magnetic tape 

Add absolute value of N in M to cleared A 

Subtract absolute value of N in M from cleared A 

Add N in M to uncleared A 

Subtract N in M from uncleared A 

Recall 50 words (one track) from drum to M 

Record 50 words in M to drum 

Add absolute value of N to uncleared A 

Subtract absolute value of N from uncleared A 

Unconditional transfer of control 

Unconditional transfer of control 

Conditional transfer of control 

Conditional transfer of control 

High speed print (Synchroprinter) 

Punch paper tape 

Multiplication (Round Off) 

Multiplication (No round off) 

Replacement of N in A to M 

Division 

Left shift 

Add associated address to number in A 

Slow speed print (Teletype) 

Place N from M in Q 

Place N from Q in M 

Drop sign of number in A 

Right shift 

Place associated address in A 

Address replacement in M 

Address replacement in M 

Half word replacement in M 

Half word replacement in M 


2 
3 
4 
5 
6 
7 
8 
9 


FD 

FE (not used) 

FF N on paper tape input to M 
800 Add N from Q to A 
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Table | 
Vocab 
| 1 AA AE 
AC 
AB AI 
| AC AE 
AD AE 
AE B: 
AF BI 
BA BC 
BB BI 
BC BI 
=) 10 BD BI 
11 BE 
12 BF Cl 
13 CA 
14 CB Cl 
15 Cl 
16 CD 
17 CE 
18 CF 
19 DA 
20 DB 
21 DC 
22 DD 
4 23 DE 
24 DF 
25 EA 
: 26 EB 
27 EC 
28 ED 
29 EE 
30 EF 
31 FA 
; 32 FB 
33 FC 
34 
35 
36 Tota 
37 
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These tapes provide the data in a convenient and flexible form to make cumula- 
tive analyses of various forms. 

Results to date include analyses of some typical Laboratory problems. In 
Table 2 the analysis of a problem in hydrodynamics is shown. Analyses are also 
given in Tables 3 and 4 for a Monte Carlo problem and for one in mathematical 
logics, respectively. A major portion of the MANIAC’s subroutines has been 
analyzed. The data for thirty-seven subroutines are shown summarized in Table 5. 
The running time for the analysis of the hydrodynamics problem was twenty-five 


TABLE 2 
Analysis of the Code for a Problem in Hydrodynamics 


Percentage Percentage 
Vocabulary Static of Static Dynamic of Dynamic 
Symbol Count Count Count Count Time 


156 13.5 3499 12.3 


S 
an 
w 
oo 


sot 


8 
0 
0 
0 
1 
93 
68 
0 
0 
0 
4 
34 
49 
13 
24 
0 
0 
93 
3 
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of Time 
AA 4 
AB 
)) AE 
BA | 
\ BB 
BC 
BD 
BE 
BF 
CA 
CB 
cc 
cD 
CE 
CF 
DA 
DB 
DC 157 1 
DD 49 
DE 28 
DF 6 : 
EA 0 
EB 118 1 
EC 104 
ED 13 
EE 11 
EF 19 
FA 38 
FB 42 
FC 0 
FD 0 
FF 17 
800 3 
Totals 1151 28333 | : 
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minutes, the Monte Carlo problem took seven minutes, and the logical problem 
took twenty minutes. Each subroutine analysis was a matter of seconds. 
Examination of these preliminary results provides quantitative estimates to 
such matters as the effect on the running time for a problem if the operation time 
of some order were decreased by some factor. In other words, attention may be 
focussed where the ‘“‘shoe pinches most.’’ For example, a decrease of fifty per cent 
in the multiplication time alone would cut the running time by 30 per cent in 
some of the typical problems. We also see that some orders could have been 


TABLE 3 
Analysis of the Code for a Monte Carlo Problem 
Percentage Percentage 
Vocabulary Static of Static Dynamic of Dynamic . Percentage 
Symbol Count Count Count Count Time of Time 
AA 37 18.2 2111 14.4 189.9 6.6 
AB 2 0.9 112 0.7 10.0 0.3 
AC 0 0.0 0 0.0 0.0 0.0 
AD 0 0.0 0 0.0 0.0 0.0 
AE 2 0.9 499 3.4 44.9 1.5 
AF 0 0.0 0 0.0 0.0 0.0 
BA 20 9.8 1326 9.1 119.3 4.1 
BB 11 5.4 730 5.0 65.7 2.2 
BC 0 0.0 0 0.0 0.0 0.0 
BD 0 0.0 0 0.0 0.0 0.0 
BE 0 0.0 0 0.0 0.0 0.0 
BF 0 0.0 0 0.0 0.0 0.0 
CA 4 1.9 89 06 . 4.4 0.1 
CB 12 5.9 395 2.7 19.7 0.6 
cc 3 1.4 231 1.5 10.3 0.3 
CD 10 4.9 786 5.3 35.3 1.2 
CE 0 0.0 0 0.0 0.0 0.0 
CF 0 0.0 0 0.0 0.0 0.0 
DA 3 1.4 909 6.2 943.5 32.9 
DB 4 1.9 356 2.4 369.5 12.9 
DC 38 18.7 2972 20.4 178.3 6.2 
DD 2 0.9 499 3.4 517.9 18.0 
DE 9 44 1085 74 146.4 5.1 
DF 1 0.4 89 0.6 6.2 0.2 
EA 0 0.4 0 0.0 0.0 0.0 
EB 9 4.4 864 5.9 56.1 1.9 
EC 7 3.4 454 3.1 29.9 1.0 
ED 8 3.9 210 1.4 8.4 0.2 
EE 4 1.9 356 2.4 67.6 2.3 
EF 6 2.9 109 0.7 8.1 0.2 
FA 4 1.9 26 0.1 2.0 0.0 
FB 0 0.4 0 0.0 0.0 0.0 
FC 1 0.4 89 0.6 7.1 0.2 
FD 2 0.9 89 0.6 7.1 0.2 
FF 0 0.0 0 0.0 0.0 0.0 
800 2 0.9 178 1.2 13.1 0.4 
Totals 201 14564 2861.8 
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sacrificed without too much pain or regret. Examination of subroutine times has 
led to a re-evaluation of some of the methods used and in some cases has resulted 
in vast improvements. 

It is planned to continue gathering statistics of problems computed on the 
MANIAC to establish a more quantitative basis for some of our ideas and 
conjectures. 


TABLE 4 
Analysis of the Code for a Problem in Mathematical Logics 
Percentage Percentage 
Vocabulary Static of Static Dynamic of Dynamic Percentage 
Symbol Count Count Count Count Time of Time 
AA 147 19.4 2025 16.5 182.2 2.1 
AB 6 0.7 20 0.1 1.8 0.0 
AC 0 0.0 0 0.0 0.0 0.0 
AD 0 0.0 0 0.0 0.0 0.0 
AE 2 0.2 4 0.0 0.3 0.0 
AF 1 0.5 1 0.0 0.0 0.0 
BA 70 9.2 1018 8.3 91.6 1.0 
BB 30 3.9 412 3.3 37.0 0.4 
BC 0 0.0 0 0.9 0.0 0.0 
BD 0 0.0 0 0.0 0.0 0.0 
BE 0 0.0 0 0.0 0.0 0.0 
BF 0 0.0 0 0.0 0.0 0.0 
CA 26 3.4 30 0.2 1.5 0.0 
CB 29 3.8 84 0.6 4.2 0.0 
co. 28 3.7 625 5.1 28.1 0.3 
CD 35 4.6 1075 8.7 48.3 0.5 
CE 1 0.1 180 1.4 453.6 5.2 
CF 0 0.0 0 0.0 0.0 0.0 
DA 0 0.0 0 0.0 0.0 0.0 
DB 1 0.1 9 0.0 9.3 0.1 
DC 90 11.9 1745 14.2 104.7 1.2 
DD 0 0.0 0 0.0 0.0 0.0 
DE 45 5.9 978 7.9 106.6 1.2 
DF 27 3.5 71 0.5 4.9 0.0 
EA 4 0.5 4 0.0 6647.6 77.1 
EB 27 3.5 1031 8.4 67.0 0.7 
EC 20 2.6 626 5.1 41.3 0.4 
ED 12 5.1 12 0.0 0.4 0.0 
EE 31 4.1 93 0.7 20.6 0.2 
EF 13 1.7 450 3.6 33.7 0.3 
FA 19 2.5 423 3.4 33.8 0.3 
FB 19 2.5 825 6.7 66.0 0.7 
FC 8 1.0 15 0.1 1.2 0.0 
FD 10 a3 ’ 30 0.2 2.4 0.0 
FE 0 0.0 0 0.0 0.0 0.0 
FF 5 0.6 15 0.1 600.0 6.9 
800 19 2.5 435 3.5 32.1 0.3 
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TABLE 5 
Analyses of the Codes for Some of the Subroutines, Summarized 
Percentage Percentage 
Vocabulary Static of Static Dynamic of Dynamic Percentage 
Symbol Count Count Count Count Time of Time 
AA 450 16.7 5659 17.6 509.3 3.5 
AB 21 0.7 123 0.3 11.0 0.0 
AC 0 0.0 0 0.0 0.0 0.0 
AD 0 0.0 0 0.0 0.0 0.0 
AE 7 0.2 23 0.0 2.0 0.0 
AF 10 0.3 10 0.0 0.9 0.0 
BA 187 6.9 2643 8.2 237.8 1.6 
BB 83 3.0 1833 5.7 164.9 1.1 
BC 0 0.0 0 0.0 0.0 0.0 
BD 0 0.0 0 0.0 0.0 0.0 
BE 0 0.0 0 0.0 0.0 0.0 
BF 3 0.1 38 0.1 3.4 0.0 
CA 66 2.4 324 1.0 16.2 0.1 
CB 68 2.5 476 1.4 23.8 0.1 
OC 92 3.4 1572 4.9 70.7 0.4 
CD 101 37 2601 8.1 117.0 0.8 
CE 42 1.5 378 1.1 952.5 6.6 
CF 0 0.0 0 0.0 0.0 0.0 
DA 15 0.5 81 0.2 84.0 0.5 
DB 15 0.5 227 0.7 235.6 1.6 
DC 349 13.0 5079 15.8 304.7 2.1 
DD 11 0.4 81 0.2 84.0 0.5 
DE 136 5.0 2184 6.8 251.1 1.7 
DF 70 2.6 429 1.3 30.0 0.2 
EA 6 0.2 6 0.0 - 9971.4 69.8 
EB 95 3.5 1634 5.1 106.2 0.7 
EC 61 2.2 1045 3.2 68.9 0.4 
ED 227 16.2 227 0.7 9.0 0.0 
EE 54 2.0 255 0.7 33.1 0.2 
EF 65 2.4 717 y 53.7 0.3 
FA 53 1.9 906 2.8 72.4 0.5 
FB 67 2.4 2214 6.9 177.1 1.2 
FC 29 1.0 36 0.1 2.8 0.0 
FD 27 1.0 46 0.1 3.6 0.0 
FF 5 0.1 7 0.0 600.0 4.2 
300 58 2.1 1147 3.5 84.8 0.5 
Totals 2572 32009 14283.3 
EvuGENE H. HERBstT 
N. METROPOLIS 
B. WELLS 
Los Alamos Scientific Laboratory 


Los Alamos, New Mexico 


Work performed under the auspices of the Atomic Energy Commission. 

1An “order” here refers to a single arithmetical or logical operation; e.g., add, multiply, 
transfer control. The ensemble of orders constitutes the computer vocabulary. An “instruction” 
refers to an order and its associated address. 

2 An interpretive routine is one which simulates a computer control and arithmetic unit. In 
effect, it seulates and performs a coded sequence of instructions. 
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Wolfram, Vega, and Thiele 

In my article ‘‘New information concerning Isaac WoLFRAM’s life and calcu- 
lations,” MTAC, v. 4, 1950, p. 185-200, special consideration is given to his 
extraordinary table of Jn x to 48 D, as published in J. C. Scnutze, Recueil de 
Table Logarithmiques, v. 1, 1778, p. 190-258. In this table are 3457 arguments, 
not 3462 as stated on p. 193, line —7, and p. 197, line 8. This erroneous statement 
was caused by overlooking the fact that in the 69 pages of the table there were 
ten groups of figures on every page except 256, where there were only nine groups. 
Thus this page reduced the estimated number of arguments by 5. Hence certain 
changes must be made in the text. 

Following the change indicated above on p. 193, for 2230, read 2225; for 928, 
read 904 ; for 533, read 552. In addition to the necessary change indicated on p. 197 
are others. First of all 9579 was a misprint for 9599. The reduction of the number 
of arguments in WOLFRAM by 5 means that there should be 79 arguments in 
THIELE, not 74, which are not in WOLFRAM; the additional arguments to the 74 
listed are : 6049, 7453, 9707, 9821, 9877. In the last seven lines of p. 196, for 3456, 
read 3457 ; for 2280, read 2225; for 74, read 79. 

In referring to VEGA’s 1794 reprint of WoLFRam’s table, p. 194-195, I failed 
to note that VEGA gave only 3451 arguments, that is, 6 less than WOLFRAM. This 
fact was brought to my attention in June 1954, by Dr. ALAN FLETCHER, of the 
University of Liverpool. I now find that the 6 WoLFRAM arguments omitted in 
VEGA are the composite numbers 2215, 2225, 2233, 2299, 2387, 2401. Since no 
one of these omissions is a prime, PETERS’ and STEtN’s Table 13, based on VEGa, 
is unchanged. 

Next, I refer to two matters in a letter of April 6, 1953, from my friend Mr. 
C. R. Cosens of the University of Cambridge. On p. 197, I had written concerning 
THIELE’s 1908 table, ‘‘Curiously enough WoLFRam’s error in no. 28 (7853) is 
corrected. This is indeed a major mystery ; the only explanation which I can offer 
is that the typesetter substituted an 8 for a 7, by mistake which THIELE did not 
observe!’’ The correction of this error in WOLFRAM was published by BURCKHARDT 
in 1817. I agree with CoseEns that a better explanation of THIELE’s achievement 
in this regard may have been that BuRCKHARDT’s correction had been brought to 
his attention. 

The second matter which Mr. CoseEns discusses at some length in his letter, 
has reference to the Caliberstabe, p. 189, lines 2-5, and footnote 34. Mr. CosENs 
shows that such a Caliper Rule, with scales, was used in connection with artillery ; 
that the weight of a round shot equals the cube of its diameter. The diameter 
(bore of the gun) would be the cube root of the weight. 


R. C. ARCHIBALD 
Brown University 
Providence, Rhode Island 


A Note on Approximating Polynomials for Trigonometric Functions 


High speed automatic digital calculators have two means available for the 
evaluation of sin x and cos x when «x is given. Either a table of values of the 
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required function may be held in the store! of the machine, and the given value 
obtained by interpolation; or the machine can calculate the required value from 
a number of terms of an infinite series. The former procedure is likely to be un- 
satisfactory if a high degree of accuracy is required since it is rarely possible to store 
enough function values to make linear, or even quadratic, interpolation feasible. 

If function values are to be calculated directly from a series it is well known 
that the ordinary Taylor series is not the best possible. Chebyshev polynomials 
have the property that » terms of their series expansion of any function define 
a polynomial of degree » which minimizes the absolute value of the difference 
between any function value and the corresponding polynomial value in the in- 
terval (—1, +1). Likewise, the Legendre polynomial expansion minimizes the 
integral square difference between function and approximating polynomial 
in (—1, +1). 

The two types of expansion find distinct applications in programming for a 
high speed computer. 


(1) If the accuracy of the work is the maximum of which the machine is 
normally capable (often 9 or 10 decimal places) the Chebyshev series is appro- 
priate, since it can give function values to the required precision with the minimum 
complexity. 

(2) If the required accuracy is less than the full capacity of the machine 
(say 3-6 decimal places) and if the results of a number of calculations are to be 
combined additively as in summing a Fourier series, then the Legendre poly- 
nomial series will be the best basis of approximation. 


In the work of this laboratory both of the above applications are of frequent 
occurrence, and since the numerical coefficients of the series required were ap- 
parently not to be found in the literature it is though€ that they may be of interest. 
to other workers in the field. 

The Chebyshev polynomials are defined, following LANczos,? by: 


To(x) = 1, Tnr(x) = cos (m x), 


and it may be shown that* the required expansions are: 


sin (xx/2) = 2 (—1)"J (4/2) T2ma1(x), 


cos (wx/2) = Jo(/2) + 2 2m (x). 
m=) 
The Legendre polynomials are conveniently defined by Ropricugs’ formula: 


Po(x) = 1, = — 1)", 
and it can be shown that +4 
sin (4x/2) = — 


cos (wx/2) = — 5+ Jsj2(4/2)- P2(x) 
+ 9+ — 
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To make use of these expansions, for finding optimum polynomials for use 
with an automatic digital calculator, it is necessary to have numerical values for 
the functions J,(2/2)Jn+1/2(#/2). Since these do not appear to have been tabu- 
lated it was thought worth constructing the table given below. Values were ob- 
tained by means of the well-known series: 


using the value ( = 2/2) taken to 20 decimal places. The resulting values were 
then rounded off to 11 decimal places and the resulting table checked by an 
application of the recursion formula: 


J,(2) = 2/2v{ J,41(2) + J,-1(2)}. 


In addition, the value of Ji:.5(x/2) was calculated directly from the recursion 
formula using the explicit value derived from the initial values: 


Jj2(/2) = 2/x, 
Jsj2(4/2) = 
and the table of values of x—* to 25 decimal places computed by GLAIsHER.* 


n Jn(x/2) n 
0 0.47200 12157 7 5 .63661 97723 7 
1 0.56682 40889 1 1.5 40528 47345 7 
2 0.24970 16291 4 2.5 .13741 70540 3 
3 0.06903 58882 9 3.5 .03212 73337 1 
4 0.01399 60398 1 4.5 .00575 32170 8 
5 0.00224 53571 2 5.5 .00083 61720 0 
6 0.00029 83476 0 6.5 .00010 23428 0 
7 0.00003 38506 4 7.5 .00001 08228 5 
8 0.00000 33522 0 8.5 .00000 10077 8 
9 0.00000 02945 7 9.5 -00000 00838 4 
10 0.00000 00232 7 10.5 .00000 00063 0 
11 0.00000 00016 7 11.5 .00000 00004 3 
12 0.00000 00001 1 12.5 .00000 00000 3 
A. D. 
University of London 
England 


and H. V. Bootu, Automatic Digital Calculators. Butterworths (London), 
1953, p. 1 
952, 7& Law Lanczos, Tables of Chebyshev Polynomials. NBS Applied Math. Ser. 9, Washington, D. C., 
1 

tEN N. Watson, A Treatise on the Theory pes Functions. 2nd ed., Cambridge, 1944, p. 21. 


BAUER, Crelle’s Journal, v. LVI, 1339, 
L. GLaIsHER, London Math. Soc., roc., 1877, p. 140. 


Continued Fraction Expansion of 2! 


The Institute for Advanced Study computer is being used to compute ex- 
tensive continued fraction expansions of certain real algebraic numbers. The 
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interest lies in comparing statistics of such expansions with known distributions 
of these statistics over random numbers. For example, KHINTCHINE! has shown 
that over random numbers x uniformly distributed between 0 and 1, the sum 
S,(x) of the first » partial quotients of x is equivalent in the sense of Bernoulli 
to Z, = n log n/log 2. 

The first result is the computation of more than 2000 partial quotients of 21. 
The table below shows S, (2+) for m = 100(100)2000, with Z, given for comparison. 
It appears that S,,(2*) oscillates considerably in relation to Z,, being most of the 
time larger, up to a factor of about 2. We do not know whether this deviation is 
significant, since! oscillations of S,(x) of this type occur for almost all x. Expan- 
sions of additional numbers, as well as more detailed statistics, will follow. 

The code depends on subroutines which do the necessary algebra on poly- 
nomials whose coefficients are p-tuples of computer words, for arbitrary and 
variable p. At present it handles cubic polynomials ; a generalization to nth degree 
polynomials is planned. The methods and results will be reported later at 
greater length. 


n n log n/log 2 S,(2!) 
100 664.4 1384 
200 1528.8 2283 
300 2468.6 2834 
400 3457.5 3471 
500 4482.9 4191 
600 5537.3 12636 
700 6615.8 18190 
800 7715.1 18777 
900 8832.4 19139 

1000 9965.8 19724 
1100 11113.6 20322 
1200 12274.6 21825 
1300 13447.6 22873 
1400 14631.7 23293 
1500 15826.1 24271 
1600 17030.2 25259 
1700 18243.2 25819 
1800 19464.8 26442 
1900 20694.4 27063 
2000 21931.6 41198 


JoHN von NEUMANN 
BRYANT TUCKERMAN 
Institute for Advanced Study 
Princeton, New Jersey 
“Metrische Kettenbruchprobleme,” Compositio Mathematica, v. 1, 1935, 
p. 361-382. 


The Values of (}) and r(3) and their Logarithms 
Accurate to 28 Decimals 


The values of I'(3), log T'(4), log were computed to 28 decimals 
using the series 


log T'(2 + x) = + — + — Cor? 
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where Cy = 1 — = lim (1 togn): Buler’s constant 


(ry = 2, 3, ---). 


= 


C, 


The values of S, = a+ytgte and y were taken from Stieltjes’ table." 


Part of the calculation was done with the assistance of Mr. E. V. HANKAM on 
an IBM (602-A type) calculating punch. Uhler’s radix table was used for getting 
the antilog of log f. The values (4) and I'(3) were required for calculating the 
power series coefficients of Bessel functions of order 4 and of functions related 
to them. 

The values were checked by the identity 


T'(3) = 2.67893 85347 07747 63365 56929 410 
T(%) = 1.35411 79394 26400 41694 52880 282 
log (3) = .98542 06469 27767 06918 71740 370 
log P(#) = .30315 02751 47523 56867 58628 174 


B. ZONDEK 
Watson Scientific Computing Laboratory 
New York, N. Y. 


1H. T. Davis, Tables of Higher Mathematical Functions. v. 11, The'Principia Press, 1935, p. 244. 


Modification of a Method for Calculating Inverse 
Trigonometric Functions 


The 605 programming that I gave recently’ fails for arguments near 2-4. The 
reason for this failure is that the double angle formulations used multiply round-off 
errors until they are intolerably large. These formulations were originally intro- 
duced to assure that cos 20 depend on both sin @ and cos 6. Upon closer examina- 
tion it was found that it is only necessary that cos 26 depend on sin 8, hence we 
may use 


cos 20 = 1 — 2 sin* @. 
The use of the above formula and 
sin 20 = 2 sin @ cos @ 
avoids the errors mentioned and is just as easily programmed for the 605. 


R. L. La Fara 


U. S. Naval Ordnance Plant 
Indianapolis 18, Indiana 


1 Ricwarp L. La Fara, MTAC, v. 8, 1954, p. 132-139. 
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REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 


1[F].—A. Table des factorisation des Nombres N* + 1 dans l’intervalle 
10000 < N < 20000. Manuscript of 84 leaves deposited in the UMT Fixe. 


This is an extension of earlier tables by the same author. Many of the entries 
have been completely factored. All unknown factors lie beyond 800 000. The 
author is preparing a table for the range 20000 < N < 30000. [For previous 
tables of this kind see MTAC, v. 2, p. 211, 252, 300; v. 3, p. 21, 118-9, 486; v. 4, 
p. 224; v. 5, p. 28, 1334; v. 6, p. 102; v. 7, p. 33-4; v. 8, p. 166.] 

T. H. SourHarp 
University of California 
Los Angeles 24, California 


2[F].—Rupotru OnprEjKA, List of the First 17 Perfect Numbers. Two typewritten 
pages deposited in the UMT FILE. 


Decimal values of the first seventeen perfect numbers, having the following re- 
spective numbers of digits: 1, 2, 3, 4, 8, 10, 12, 19, 37, 54, 65, 77, 314, 366, 770, 
1327, 1373. 

Computation of the first twelve perfect numbers was done by the author with 
the use of his table of 2", m = 1(2)411. The last five were computed by H. S. 
UHLER and were checked by the present author. The present list is believed by 
the author to be error-free. 

T. H. SourHarD 
University of California 
Los Angeles 24, California 


3[F].—R. J. Porter, A List of Groups and Series to serve for computations of 
Irregular Negative Determinants of Exponent 3n. 274 typewritten pages de- 
posited in the UMT Fie. 


This is very closely related to the author’s UMT 155 [MTAC, v. 7, p. 34] 
and UMT 185 [MTAC, v. 8, p. 96-7]. 


T. H. SouTHAaRD 
University of California 
Los Angeles 24, California 


4[K].—D. V. LinpLtey & J. C. P. MrterR, Cambridge Elementary Statistical 
Tables. Cambridge University Press, London, 1953, 35 p., 21.9 X 27.9 cm. 
Price $1.00. (paper) 


“This set of tables is concerned only with the commoner and more familiar 
and elementary of the many statistical functions and tests of significance now 
available-—The more familiar statistical tests are either based directly on the 
normal distribution or, in the case of the #, x? and F tests, they are derived there- 
from. Percentage points for these tests are provided in the tables, mainly for 
significance levels 5%, 1%, and 0.1% in both one-sided and two-sided tests.— 
Tables of the more common transformations (of the data), square root, logarithm, 
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inverse circular and hyperbolic root-sines, together with that for the correlation 
coefficient, have been included.”” (From the Preface.) 
A list of the tables follows: 


1, 2. Cumulative normal (x) to 4D with differences for x = 0(.01)3(.1)4 
with normal ordinates to 4D for x = 0(.1)4, special values of x for (x) 
= .001 (.001).03(.002).05(.05).50 and several special values. 

3. Percentiles 87.5, 95, 97.5, 99, 99.5, 99.9, 99.95 of the #distribution for 
degrees of freedom 1(1)10, 12, 15, 24, 30, 40, 60, 120, © to 2D. 

4. z = tanh’? for r = 0(.02).8(.01).94(.001)1, to 3D with first differences. 
This is a transformation of the correlation coefficient. 

5. Percentiles .5, 1, 2.5, 5, 90, 95, 97.5, 99, 99.5, 99.9 of the x? distribution for 
degrees of freedom 1(1)30(10)100 to 3 or 4S. 

6. Conversion of range to standard deviation for sample size n = 2(1)13 to 
4D. (Ratio of expected values.) 

7. Percentiles 95, 97.5, 99, 99.9 of the F-distribution for degrees of freedom 
v; = 1(1)8, 10, 12, 24, 2 and degrees of freedom vz = 1(1)30(2)40, 60, 120, « to 
3 or 4S (in general 2D). 

8. 2000 random digits. 

9. Various functions. 

For n = 0(1)100: n?, Vn to 4D, 1/n to SD, 1/Vn to SD. 

For x = 0(.01)1: sin Vx, sinh Vx, sinh V10x, sinh- ¥100x all to 3D with 
first differences. ee 

For x = 0(.01)10: x* exact; Vx, ¥10x, 1/x, 1/Vx, 1/Vi0x each to 4S when the 
first significant digit is >2 and to 5 figures when <2, with first differences; log x 
to 4D with first differences. 

For log ¢ = 0(.001)1:¢ to 4D for ¢ < 2 and to 3D for ¢ > 2 with first differences. 

10. log 2! for m = 0(1)300 to 4D. 


The tables are arranged in a convenient format with notes on interpolation 
and asymptotic expressions for values beyond the given tables. A simple but not 
exact description of the tables is that this book is an abbreviated form of the 
FIsHER & YaTEs Tables,! since that book contains, among other tables, most of 
those listed above. The principal area in which the book under review is more 
complete is in Table 7 and part of Table 9. Tables 3, 5, and 7 are based on tables 
from Biometrika but contain some additional values. Values which are reported 
here that do not appear either in the Biometrika tables or those of FisHER & 
Yates are those for vz = 32(2)38 for all percentiles and those for »; = 7 and 10 
for percentile 99.9. A number of differences of a single unit in the final place were 
noted in the 99.9% F-table between the FisHER & YATES Tables! and the tables 
under review. 

W. J. Drxon 


University of Oregon, 
Eugene, Ore. 


1R. A. Fisher & Frank YATES, Statistical Tables for Biological, Agricultural, and Medical 
Research. London & Edinburgh, 4th ed., 1953. 
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5[K].—H. G. Romi, 50-100 Binomial Tables. New York, John Wiley & Sons, 
1953, xxvii + 172 p., 18.7 XK 22.4 cm., $4.00. 


These useful tables supplement the National Bureau of Standards Tables! 
which give individual terms and partial sums of terms of (g + p)" to 6D for 
b = .01(.01)50 and m = 2(1)49. Here the same quantities are given also to 6D 
for m = 50(5)100. The arrangement differs from that of the previous tables by 
giving all the entries for each n, p pair in adjacent columns, values of individual 
terms in one and cumulative sums in the other. 

Dr. RomiG provides very adequate explanation and illustrative examples to 
assure proper use of the Tables. The author also includes exact interpolation 
forniulas determined by proper transformations of the interpolation formulas for 
the incomplete beta-functions. These formulas may be used for those probability 
determinations for intermediate p and m values for which ordinary linear inter- 
polation is not satisfactory. In addition, there is a satisfying list of references 
which more completely cover the theory of interpolation. 

These tables will find many uses in numerous areas of statistical analysis. 
Especially will they be useful in the area of statistical quality control where more 
refined evaluations of the binomial probabilities are necessary to determine the 
protection afforded by proposed sampling procedures in process control tech- 
niques as well as in acceptance plans. 


G. W. McELRAtTH 


University of Minnesota 
Minneapolis, Minnesota 


1NBSCL, Tables of the Binomial Probability Distribution. AMS no. 6, Washington, 1950. 
[MTAC, v. 4, p. 208-209.] 


6([K].—W. L. Stevens, “Tables of the angular transformation,’ Biometrika, 
v. 40, 1953, p. 70-73. 


The following form of the angular transformation 


@ = 50 — Xarcsin (1 — 2p), » = V1000 


has been tabulated by the author in order to provide ‘‘a table similar in accuracy 
to that of the table of probits given by FrsHER and Yares.’’! The author suggests 
that the present form of the transformation has the following advantages: (i) 6 
ranges from 0.327 to 99.673 and therefore has almost the maximum possible 
accuracy for any given number of significant figures; (ii) the weight is given by 
the extremely simple expression /1000, where n is the number of observations; 
(iii) like the probit function, complementary values of the function correspond to 
complementary values of the argument p = 50% giving @ = 50. 

Three tables are presented, each to 3D. Table 1 gives @ for 100p = 0(0.1)50 
while for 100p > 50 entry is made for 100(1 — ) and the transformed value is 
equal to 100 minus the tabular value. Proportional parts are given for linear 
interpolation. For small percentages Table 2 gives @ for 100p = 0(0.01)2.0, with 
proportional parts for linear interpolation when .05 < 100 < .2. For 100p < 98, 
6 is determined as above by subtracting the tabular value obtained in Table 2 
for 100(1 — p) from 100. For 100p < .05 the formula @ = 0.327 + 6.325Vp may 
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be used. Table 3 provides the values for proper fractions whose denominations 
are less than or equal to 30. 


S. B. Lirraver 
Columbia University 
New York, N. Y. 


1R. A. FisHEeR & F. Yates, Statistical Tables for Biological, Agricultural and Medical Research. 
London & Edinburgh, 1938. 


7[K ].—A. C. Conen, Jr. & Jonn Woopwarp, “Tables of Pearson-Lee-Fisher 
functions of singly truncated normal distributions,” Biometrics, v. 9, 1953, 
p. 489-497. 


A normal variable, x, with frequency function, 
o(x) = exp (x — m)*/2c?, 
for all real x, truncated (below) at xo’, gives the truncated variable, x’, with fre- 


quency function, f(x’) = $(x’)/Io(x’), xo’ < x’, where Io(x’) = f . o(x)dx. It is 


required to estimate m and o on the basis of a sample of observations on 
x =x! — xq. 

Pearson and Lee! gave estimates based on the first two observed moments. 
FisHER’ showed their estimates to be maximum likelihood estimates. The authors 


estimate § = — and oa by the relations (equivalent to the earlier ones) 


and 


- = ne, 


where z = $(£)/Io(é). To facilitate computations they give tables of A(£) and 
g(t) to 8D except for the largest values of § (where 7D and 6D are given) for 
§ = — 4.(.1) —2.5(.01).5(.1)3. The authors suggest using (1) to estimate = and 
(2) with = from (1) to estimate ¢. 

The variances of the estimates and the correlation coefficient between the 
estimates are given by 


1 — 2(z — &) 


2 — &(z — &) w'(é) 


var(¢) = 


W'(é), 


var(¢) = 


s-—€ 


p(é, ¢) = 
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Tables of W’(é) and w’(£) are given to 6D and, of p(é, «) to 4D for = —3.(.1)2, 


Katz 
Michigan State College 
East Lansing, Michigan 


Kari Pearson and Atice LEE, the generalized probable error in multiple normal 
Biometrika, v. 6, 1908, 
?R. A. FisHer in BAAS Math. ables, v. I, London, 1931, p. xxvi-xxxv. 


8[K].—P. V. K. Iver & A. S. P. Rao, “Theory of the probability distribution of 
runs in a sequence of observations,’ Indian Soc. Agricultural Stat., Jn., v. 5, 
1953, p. 29-77. 


This paper investigates the theory of runs in which the succession of observa- 
tions may be equal, ascending or descending. As such it augments the existing 
work on runs in which the equality assumption was not allowed. The purpose of 
the paper is to investigate the distribution of the number of ascending, descending 
and stationary runs in a sequence of m observations. Both the infinite case, where 
a particular value has a given probability of occurrence, and the finite case, where 
one knows the number of times a given value has occurred, is considered. For each 
of the three types of runs, the various related configurations are tabulated along 
with their probability of occurrence in the sequence and the number possible. 
No actual distributions are attempted in the paper. Variances and covariances 
for k kinds of elements with equal probabilities of occurrence are tabulated to 4D 
for K = 2(1)5, 10 and n = 30, 40, 50, 75, 100 for ascending, stationary, and 
descending runs as well as for the total number of runs. The author also lists the 
algebraic expressions for the covariances of runs of lengths p and g for p = 1(1)4, 
q = p(1)5, except for p = 4 only g = 4 is given; of runs of length » and q or 
more for p, g = 1(1)4; and of runs of lengths p or more and g or more for p = 1, 
2, 3, g = p + 1(1)4. For junctions the actual distributions are given for values 
of m = 4(1)?7. 

Mention is made of the possibility of using the results for testing randomness 
but the actual discussion of such applications are to be given in a separate article. 

Cart F. Kossack 
Purdue University 
Lafayette, Ind. 


9[K].—S. H. Aspet-Ary, “Tables of generalized k-statistics,” Biometrika, v. 41, 
1954, p. 253-260. 


The author gives a complete table of order 12 for expressing the sample &,.... 
statistics of TUKEY! in terms of the augmented monomial symmetric functions 
of Davip & KENDALL? and vice versa. This is a very considerable extension of 
the table of WisHART,* which were through order 6, since results for lower orders 
are at once obtainable by a simple rule for the deletion of subscripts. 

Gc..c. 


. W. TuxeEy, “Some sampling simplified,” Amer. Stat. Assn., Jn., v. 45, 1950, p. 501-519. 
. W. Davin & M. G. vey ~ Tables of eymmetric functions—part I, Biometrika, v. 36, 
1949, 431-449. [MTAC, v. 4, Pon. 
WisHaRT, “Moment coefficients the k-statistics in samples from a finite population,” 
Biometrika, v. 39, 1952, p. 1-13. [MTAC, v. 7, p. 97.] 
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10[K ].—J. H. Capwe “The statistical treatment of mean deviation,” Bio- 
metrika, v. 41, 1954, p. 12-18. 


The author wants to obtain properties of the mean deviation m, which are 
analogous to the well-known properties of the standard deviation o of a normal 
distribution. To this end, the distribution of the quotient m/o is approximated by 
the x? distribution. He matches the first two moments of the two distributions 
with a small discrepancy for the third moment. Let m(k, m) be the average of k 
mean deviations, each for samples of size from a normal population with stand- 
ard deviation ¢. Then it is shown that c[m(k, n)/o }'* has approximately the x* 
distribution with v degrees of freedom. Table 1 gives the values c to 4S and » to 
1D, the expected values of m(k, n)/o to 4D and the variances of m(1, m)/o to 5D 
as functions of k and m for k = 1(1)i0, m = 4(1)10. Table 2 gives the same values 
for k = 1(1)5, m = 10(5)50. Table 3 gives the lower and upper 2.5 percent and 
5 percent points of the probability function of m(k,m) for k = 1(1)10 and 
n = 4(1)10 to 3D. The values are exact for k = 1, and for other values the error 
will not exceed .003. For values of k beyond 10 a normal approximation can 
be used. 


E. J. GUMBEL 
Columbia University 
New York, N. Y. 


11[K].—H. Weer, “A new type of control chart limits for means, ranges, and 
sequential runs,”’ Amer. Stat. Assn., Jn., v. 49, 1954, p. 298-314. 


In the usual theory of quality control charts the probability of a sample point 
falling outside the control limits, when the system is in control (a Type I error) 
is kept constant regardless of the number of items involved at each sample point. 
In the author’s theory the average number of false alarms (Type I error) is a 
fixed percentage of the number of items tested, independent of the sample size. 
Given the variable X, normally distributed with known mean, m, and standard 
deviation, a, let p be the probability that a random sample of m causes a false 
alarm at the upper limit, let a be the average number of articles tested before the 
false alarm is raised; then a = n/p. In Table I the author lists the values of B 
and B/Wn to 2D, for a = 5000, n = 3(1)6, 8, 10, since the upper control limit is 
given by m + Bo/Vn. The lower control limit is given by m — Bo/¥n. In Table II 
are listed the values of B/Vn to 3S for n = 3(1)10(5)50 and a = 1000(1000)5000. 
Suppose the population mean m changes to m + ko, k > 0. For a given n and B, 
the average number of items tested A (m) is a function of k. A(m) is plotted against 
k for a = 2000, n = 5, 10, 20, 50 (Chart I), and for a = 5000, same values of n 
(Chart II). From these charts the author concludes that if large values of k are 
expected, k > 1.6 say, then small samples, e.g., m = 5, should be used; while if 
small values of k are expected, k < 1 say, large samples of m, e.g., m = 10 or 20, 
are more economical. In Chart III A(m) is plotted against & for fixed sample size 

= 10, a = 1000, 2000, 3000, 5000. This chart indicates that small values of a 
are useful only for the detection of small changes of the population mean. 

Similarly for the control chart of the range in Table III are given to 2D 
the values of Wio < R < Woe, the control limits for the sample range R, for 
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n = 3(1)10, a = 1000(1000)5000. In Chart IV A(m) is graphed for n = 5, 10, 
1.0 < k < 2.4, for the range. In Chart V A(m) is graphed for nm = 10, a = 1000, 
2000, 5000 for the range. For the range it is assumed that o changes from ¢ to 
= Ko, K>1. 

The author next turns to the use of runs for controlling the mean, where 
X = length of run. In Chart VI is graphed the value of A(n) for 0 < k < 1.4, 
n = 4, 10, 20 and a = 4000, \ = 2, and in Chart VII A(m) is graphed for the 
same values of m and a but A = 3. In Chart VIII A(m) is graphed 0 < k < 1.3 
for \ = 1, m = 20, and for \ = 3, m = 4. In the range 0 < k < 1, the use of 
\ = 1, m = 20 is superior from the power sense to A = 3, n = 4. 

The use of runs for the range charts reduces rather than improves the power 
of the chart and hence no charts are given. Considerable use is made of the author’s 
two previous papers.'? 

L. A. AROIAN 
Hughes Aircraft Company 
Culver City, California 


2H. WEILER, “The use of runs to control the mean in quality control,” Amer. Stat. Assn., 
Jn., v. 48, 1953, p. 816-825. 


12[K].—S. RosensauM, ‘Tables for a nonparametric test of location,’ Annals 
Math. Stat., v. 25, 1954, p. 146-150. 


If two independent random samples of sizes m and m are drawn from a con- 
tinuous statistical population, the probability that exactly 7 points of the sample 
of m will exceed the greatest value of the sample of ” is: Q; = n(7)B(n + m — i, 
4 + 1), where B is the complete Beta function. We can fix a probability level « 
and determine a value so such that: 


oF e< Qi. 
#=0 #=0 


This paper presents tables of s = so +1 for « = .99,.95 and m, m = 1(1)50. 
These results can be used to test whether two samples came from the same popu- 
lation. The argument is identical if the number of values of the sample of m 
which are less than the smallest value of the sample of n is considered. 

J. E. 
Lockheed Aircraft Corporation 
Burbank, California 


13[K ].—H. RusBen, the moments of order statistics in samples from normal 
populations,” Biometrika, v. 41, 1954, p. 200-227. 


Many statisticians have treated the problem of evaluating the moments of 
order statistics in small samples drawn from normal populations. The results 
have been fragmentary, partly because no one has developed a systematic ap- 
proach to the problem. The author gives a systematic approach in this paper. The 
method involves showing from a geoinetrical point of view that the moments of 
normal order statistics as well as the moment generating function of any order 
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statistic are closely related to the volumes of members of a class of hyperspherical 
simplices. These volumes involve calculation of a function ag(x) (see the paper for 
details). Table 1 gives values of this function for x = 2(1)12 and B = 1(1)49 to 
8 to 10D. These values are used to compute Table 2, which gives the first ten mo- 
ments of the extreme members (smallest or largest values) in samples of size n, 
1 <2” < 50 to 9 or 10S. Table 3 gives the second, third, and fourth moments 
about the mean of extremes in samples of size 1 < n < 50, as well as the standard 
deviations, together with 6; and 6, — 3 to 7 or 8D. 
BENJAMIN EPSTEIN 

Wayne University 

Detroit, Michigan 


14[K ].—T. J. Terpstra, ‘The exact probability distribution of the T statistic 
for testing against trend and its normal approximation,’ K. Ned. Akad. van 
Wetensch, Proc., s. A, v. 56, 1953, p. 433-437. 


Let xa, i = 1,---,1,4 = 1, ---, be random samples of the random vari- 
ables x; which the null hypothesis states to have identical frequency distributions. 
Let Witcoxon’s statistic,! be the number of pairs (h, k), h < ni, k < ny, 
with i < jand xi, =1,---,l,h,k =1,---, m3 then T = Die; Vu, 
a generalization of the T defined by Mann & WHITNEY,’ was studied by the 
author® as a test against the alternate hypothesis of an upward term in the x;. 
In the present paper, in Table 1 the author gives the exact distribution of T to 
3D for m; < mz < m3 < 5, and in Table 2 gives the .005, .01, .025, .05 and 0.1 
significance levels for T (the smallest value of T for which the probability of its 
being exceeded is no greater than the relevant significance level) for the same 
values of the n’s. The values are also given for the normal approximation which 
are such that they indicate that for n; > 5 this approximation is good. 

1F. Witcoxon, “Individual comparisons by ranking methods,” Biometrics Bull., v. 1, 1945, 
FUT B. Mann & D.R. Wuirney, “On a test of whether one of two random variables is sto- 
chastically larger than the other,” Ann. Math. Stat., v. 18, 1947, p. 50-60. 

*T. J. Terpstra, “The asym ch arid and consistency Kendall's test against trend, 


when ties are present in one K. . Akad. van Wetensch., Proc. s. A, v. 55, 1952, 
p. 327-333. 


15[K ].—H. O. Hartiey & H. A. Davin, “Universal bounds for mean range and 
extreme observation,” Annals Math. Stat., v. 25, 1954, p. 85-99. 


The authors extend the theory of universal upper and lower bounds for E(w,), 
and universal upper bounds for E(x,), where x, is the standardized extreme 
variate and w, the standardized range of a sample of m. Table I gives the upper 
bound of E(x,) to 4D for any ,population for m = 2(1)20. Also given for com- 
parison are previously known values for symmetric populations. Table II gives 
the universal lower bound for E(w,) over distributions with finite range 
—X<x< X. For m = 2(2)12, the bound is given to 3D for X = 1(1)5; for 
n = 12(2)20, it is given to 3D for p = .95(.01).99, where p = X?/(1 + X?). 
Specific values are easily computed from equation (58). As X — «, E(w,) — 0 
in agreement with previous results. It is shown that universal upper bounds 
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previously computed for the case of symmetric populations with infinite range 
are also applicable to the non-symmetric case, and to the case of finite range 
unless X < X,, where X, is the limit of the range for the population which 
maximizes E(w,). Values of X, to 3D are tabulated for m = 2(1)20, and the 
algebraic form for the bound when X < X, is given. 


C. A. BENNETT 
General Electric Company 
Richland, Wash. 


16(K].—R. E. Becuuorer, “A single-sample multiple decision procedure for 
ranking means of normal populations with known variances,” Annals Math. 
Stat., v. 25, 1954, p. 16-39. 


Let Xi, ---,X, be independent normal random variables of unit variance. 
Table I gives to 4D the value of d such that y = Pr{X., ---, Xen < Xe 
+d,---,X;%+d} for y = .05(.05).8(.02).9(.01).99, .999, .9995 and for k = 2(1)10, 
t = 1(1)[R/2] as well as 10 other pairs (k, #). Table II gives to 4D the value of d 
such that y = Pr{X; < X2+d<X3+ 2d} for y = .2(.05).8(.02).9(.01).99. 
The tables enable one to compute the numbers of observations needed from nor- 
mal populations of known variance in order to have confidence y in certain 
statements about the order of the population means. 

J. L. Hopcegs, Jr. 


University of California 
Berkeley, California 


17[K].—R. E. BEcHHOFER & MILTON SoBEL, “A single-sample multiple decision 
procedure for ranking variances of normal populations,’ Annals Math. Stat., 
v. 25, 1954, p. 273-289. 


Let U, V, W, X be independent chi-square random variables, each with n de- 
grees of freedom. The paper provides 5D tables of Pr{ U < @V}, Pr{U < @V, 6W}, 
P{U, V < 0W}, Pr{U < 6V < @W} and Pr{U < @V, @W, 0X}, for m = 1(1)20 
and @ = 1.2(.2)2.2. The tables provide the confidence coefficients of certain 
statements about the order of the variances of normal populations. 

J. L. Hopcgs, Jr. 
University of California 
Berkeley, California 


18[K].—C. W. Dunnett & Mitton Soset, “A bivariate generalization of 
Student’s ¢-distribution with tables for certain special cases," Biometrika, 
v. 41, 1954, p. 153-169. 


The authors consider the simultaneous distribution of two variates, t; = 2;/s 
and t; = 22/s. The z; follow a normal bivariate distribution with zero means, the 
same variance o”, and correlation p. The variance o? is assumed independently 
estimated by s* with m degrees of freedom. The probability integral is: 


Prob {ty < h; te h} = P, 


Tables of P and h are presented for nm = 1(1)30(3)60(15)120, 150, 300, 600, « 
and p = .5and —.5. P is given to 5D for h = 0(.25)2.50 and 3.00, plus some addi- 
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tional values for larger h when 2 is small. hk is given to 3D for P = .50, .75, .90, 
.95 and .99. 
An asymptotic expansion is derived for P and h; 


4 
P=>Ai/n* and h= > 
i=0 i=0 


Values of the A; and B; to 6D are presented for the same values of h and P 
mentioned above. 
This distribution has applications in certain multiple decision problems. 


R. L. ANDERSON 
North Carolina State College 
Raleigh, N. C. 


19[K ].—J. M. Sencupta, “Significance level of }°x*/(0x)* based on Student’s 
distribution,” Sankhyd, v. 12, 1953, p. 363. 


In testing whether the mean of a normal population sampled is equal to a 
given value y, one ordinarily applies Student’s “‘t,” i.e., t = ¥n(g — u)/s, where 
N N 

1 

(yi — is the usual unbiased estimate of popu- 

lation variance. The author of the paper reviewed has noted that if one puts 
x; = Yi — p, i.e., considers deviations about the hypothetical mean tested, then 


4a) _f+N-1 
(2x)? Ne 


and hence percentage points for the right-hand member of (1) can be obtained 
from percentage or significance levels of Student’s ‘‘t.”" A table of such percentage 
points would therefore simplify computations for the test of a hypothetical mean 
and will in fact be very useful and time-saving when many f-tests are to be carried 
out. The presented table contains the 5% and 1% significance levels to 4D for 
F+N-—1 
oo for sample sizes N = 2(1)30(10)60, ©. The user of the table should 
note that significant values of (# + N — 1)/Nf are less than those tabulated. 
F. E. Grusss 


I= wands = 


Ballistic Research Laboratories 
Aberdeen Proving Grounds, Md. 


20[K ].—J. C. Sprrz, “Matching in psychology,” (Dutch, English Summary), 
Statistica, v. 7, 1953, p. 23-40. 


Let p, be the probability of exactly r matches in a random matching of two 
similar decks of 3 distinct cards. Then, po = 4, p1 = 4, b2 = 0, ps = 3. Ina 
series of nm such random matchings, let R = r; +---+ 1, be the total number of 
matches (0 < R < 3m) and let P,,, = Pr(R > a). Obviously, 


(1) Pas + + 


from which the value P,,, is tabulated to 3D for m = 1(1)30 and a = 0(1)3m. 
If in an actual experiment the number R = a of matches is such that P,,. < .05 
(say), one rejects the hypothesis that the matchings were random. 
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For large n, by the central limit theorem, R is approximately normal. Now, 
R has a mean p» = 2, a variance o? = m and a skewness y; = 1/ Vn. Thus, for 
large n,t = (R — 4 — n)/Vnis about N(0, 1). For m = 30, the resulting approxi- 
mation to P,,. appears to be fairly good. An even better approximation is ob- 
tained by replacing ¢ by a type III variable with » = o = 1, y; = 1/Vn; then 
P,,, can be readily computed from Sa.vosa’s tables.! 

The reviewer would expect a good approximation to P,,, by replacing R by 
a Poisson variable with parameter » which has up = o? = n, y; = 1/¥Vn. 


J. H. B. KEMPERMAN 
Purdue University 
Lafayette, Ind. 


ee = Satvosa, “Tables of Pearson’s type III function,” Ann. Math. Stat., v. 1, 1930, follow- 
ing p. 


21[K].—J. W. distribution of the difference in total rank value 
for two particular objects in m rankings of m objects,’ British Jn. of Stat. 
Psychology, v. 7, 1954, p. 45-49. 


Let riz (¢ = 1, ---,m; 7 = 1, be m rankings of m individuals, where 
the total ranks assigned to two particular individuals a and 6 are of interest. 
The author considers the statistic d = ¥ (ria — re), the difference in total rank 

i=1 
values, and obtains the one-sided cumulative distribution function $P[|d| > k] 
on the assumption of randomness to 5D for n = 2, 3, n = 3(1)8;m = 4, n = 3(1)7; 
m = 5,n = 3(1)5;m = 6, = 3,4; m = 7, 8, n = 3. The first four moments of 
d and a normal approximation are also given. 2 


INGRAM OLKIN 
Michigan State College 
East Lansing, Mich. 


22[K].—Ws. R. THompson, Tables of the Four Variable N- and Psi-Functions. 
9 + 88 typewritten pages (ozalided), deposited in the UMT FILE. 


Letn =r+s,n’ =r’ +s’ and define 


N(r,s,7', 8’) = s, 7’, s’) 


tss,r’ , , 1 
t=0 


where is the binomial coefficient. 


These tables give exact values of N for appropriate positive integral arguments 
such that 3 < m < 20, 2 < mn’ < 20, nm > n’ and 7S approximations (in some 
cases exact) of 10’y for appropriate non-negative integral arguments such that 
1<”<20,0< <20,”>n'’. 

The word appropriate here refers to the fact that a user of the table may 
take advantage of identities which result from permuting the arguments in certain 
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ways. A nine-page explanation of the tables and their method of construction 
is included. 


The tables include, as subheadings in each ‘“‘cell’’ of the table, parenthesized 
). The author believes the tables to be 


vaiues of n, n’, and D ie +2 
n+1 
error-free. 
The author’s address is: Division of Laboratories and Research, N. Y. State 
Department of Health, Albany 1, New York. 


T. H. SourHarD 
University of California 


Los Angeles 24, California 


23[L].—NBS Applied Mathematics Series, No. 32, Table of Sine and Cosine 
Integrals for arguments from 10 to 100. U. S. Government Printing Office, 
Washington, D. C., 1954, xvi + 187 p., 20.5 X 27 cm. Price $2.25. 


The first edition! of this volume appeared in 1942. In the present second edition 
the bibliography has been brought up to date, the table of (1 — ») has been 
replaced by a table of p(1 — »)/2, the table of E, and F, has been added. 

The principal table (180 p.) gives 10D values oi 


Si(x) = fa, = as 
0 t t 

with second central differences for x = 10(.01)100. Interpolation by EvERETT’s 
formula gives accuracy to within 1.2 units of the tenth decimal place. 

Auxiliary tables: mx/2 for n = 1(1)100; 15D. p(1 — )/2 for p = 0(.001)1; 
exact values. The EvEREtTT coefficients E,(p) and F2(p) for p = 0(.001)1; 7D. 

In the Introduction A. N. Lowan gives the fundamental formulas for these 
functions, the method of computation and the preparation and checking of the 
manuscript; and describes both direct and inverse interpolation and their accu- 
racy. A bibliography of tables, references to applications, and a list of texts and 
handbooks is also given. 

There are graphs of Si(x) and Ci(x), a Preface by A. V. Astrn, and a Foreword 
by J. A. StRATTON. 

A. E. 
1NYMTP, “Table of sine and cosine integrals for arguments from 10 to 100." New York, 1942. 


24[L].—S. Rusuton, ‘On the confluent hypergeometric function M(a, 7, x),” 
Sankhya, v. 13, 1954, p. 369-376. 

S. Rusuton & E. D. Lane, ‘Tables of the confluent hypergeometric function.” 
Sankhyd, v. 13, 1954, p. 377-411. 


The second of these papers gives 7S values of 


— + j) x? 

(18) = 
for y = .5(.5)3.5, 4.5; x = .02(.02).1(.1)1(1)10(10)50, 100, 200, and an extensive 
range of integer and half-integer values of a. This range varies between 0 < a < 25 
and 0 < a < 50, and the interval is .5 or 1. 
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In the first paper the author describes the principal properties of the confluent 
hypergeometric function, its application to facilitate certain sequential tests of 
composite hypotheses, and the construction of the tables. When a = y or y + 1, 
the confluent hypergeometric function can be expressed in terms of the exponen- 
tial function, and from here the recurrence relations enable the computer to 
proceed to other values of a, as long as a — vy is an integer. In other cases, the 
power series expansion was used for x < 5, and the asymptotic expansion, im- 
proved by Airey’s converging factor, for x > 5. 

Companion tables' for y = 3, 4 were reviewed in RMT 1003 (MTAC, v. 6, 
1952, p. 155-156). 

A. E. 

1P, Natu, “Confluent hypergeometric function,” Sankhyd, v. 11, 1951, p. 153-166. 


25[L].—Mitton Asramowi1z & Puitip ‘Evaluation of Coulomb 
wave functions along the transition line,” Phys. Rev. (2), v. 96, 1954, p. 77-79. 


Work on Coulomb wave functions at the NBS Computation Laboratory has 
been reviewed in RMT 1091 [MTAC, v. 7, p. 101-102], UMT 186 [MTAC, v. 8, 
p. 97], and RMT 1249 [MTAC, v. 8, p. 224]. In continuation of this work, the 
authors derive asymptotic expansions of Fo, Fo’, Go, Go’ for p = 2n, in descending 
powers of 8 = (%n)'. They also indicate the computation of Fy and Gp» for p near 
2n by means of Taylor series. 

The paper contains 7D tables of Fo, Fo’, Go, Go’ for p = 2n = 0(.5)20(2)50. 
The tabular values were computed to 9D on the SEAC of the NBS by means of 
a program prepared by C. E. Fréberg and based on numerical evaluation of 
integral representations. The 7D given in the tables are stated to be correct to 
within one unit of the last place, and five-point Lagrangian interpolation will 
yield full accuracy for p > 3. 

In a companion paper,' expansions of Coulomb wave functions (for any L) 
are obtained for the region 0 < p < 2p; = 2n + 2[y? + L(L + 1) ]. 

A. E. 


1 Mitton Aspramowitz & H. A. ANnTosiIEwicz, “Coulomb wave functions in the transition 
region,” Phys. Rev. (2), v. 96, 1954, p. 75-77. 


26[L].—J. Ciunte, “On Bose-Einstein functions,’ Phys. Soc. Proc. Sect. A, 
v. 67, 1954, p. 632-636. 
The author provides formulas useful for the computation of the function 
dz 
(the Cauchy Principal Value of the integral is to be taken when > 0): these 


formulas include asymptotic expansions for large || and power series expansions 
for small ||, and the formula 


(1) Gi(n) = 


(2) Gr(n) = + (2"n) 


where F;, is the Fermi-Dirac function.' 
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4D numerical tables of G;(m) are given for 7 = — 3(.2) — .6(.1).6(.2)20: out- 
side this range the asymptotic formulas are valid. The values were computed 
from (2) and available tables' of F;, they were checked by differencing and, in 
the interval —.5(.1).5, also by independent computation from the power series; 
and it is stated that the error in the present tables is less than 1 unit in the fourth 
decimal place. 

The power series expansion of G; and graphs of multiples of G_;, G;, Gy were 
given by ROBINSON.” 

A. E. 


1 McDouecatt & E. C. STONER, Soc. Phil. Trans., v. 237A, 1938, p. 67-104. 
2]. E. Ropinson, Phys. Rev., s. 2, v. 83, 1951, p. 678-679. 


27[L]}.—James G. Berry, Tables of Some Functions Related to the Legendre Func- 
tions P,,-™(x) and Q,(x) when n is a complex number. Two copies, each 28 pages 
of copied typescript, deposited in the UMT Fie. 


—m/2 d 
) (x)], for 


P,-"(x) and of ( 


Tables of ( 


m = 0(1)20, and of Q,(x) and Ps [Q,(x)], where P,,-"(x) is the associated Le- 


gendre function of the first kind and Q,(x) is the Legendre function of the 
second kind, for m = —0.5 + 10.24595735 + i(10.18477501). The range of x is 
.000 (.032).960(.002).998, and x = .999875. 

Real and imaginary parts of the two functions are in most cases given to 8S, 
although only spot checks were made and the author admits some instances in 
which the error is +5 in 7th S. Computation was done by a modification of the 
Runge-Kutta method on MIDAC, in connection with the author’s Ph.D. Disser- 
tation in Engineering Mechanics. 

T. H. SouTHARD 


University of California 
Los Angeles 24, California 


28[L ].—J. Bercuuts, A Table of Some Integrals. R 245, Mathematisch Centrum, 
Amsterdam. Eight page mimeographed typescript, deposited in the UMT FILE. 


This table contains 


8D values of f,(x) = tanvdv, mn =41(1)5, x = 0(0.05)1.50 
0 

8D values of g,(x) = cot ode, m =1(1)5, x = 0(0.05)2.50 
0 

7D values of F,(x) = tanh odo, m = 1(1)4, x = 0(0.02)1.98 
0 


7D values of G,(x) = f *“cothvdy, m=1(1)4, x = 0(0.02)1.98. 


E. 
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these | 
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The error is stated to be 10-* in the last 3 functions, 3 K 10-* in f,(x). F(x) 
and G,,(x) were computed by the ARRA;; the other functions were hand computed 
and checked by differencing on a National Accounting machine class 31. 


29[L].—Staff of the Computation Department of Mathematisch Centrum, 
Amsterdam, Table of Polylogarithms, Report R24, Part I: Numerical Values. 
53 mimeographed pages deposited in the UMT Fine. 


In three tables are given 10D values of F,,(z), defined by ~ h-*2* for |z| <1 


and for other values of z by analytic continuation, called pelylogarithos of order 
nm and argument z. In all 3 tables, m = 1(1)12. 


Table I: z2=x;x = —1(0.01)1 
Table Il: 2 = ix; x = 0(0.01)1 
Table III: 2 = e*#/?; q@ = 0(0.01)2. 


The maximum error is stated to be 10-"°. 


30[S].—Ann T. News, “Graphs of the Compton energy-angle relationship and 
the Klein-Nishina formula from 10 Kev to 500 Mev.,” National Bureau of 
Standards Circular 542, 1953, iv + 89 p. 


This Circular contains eight principal graphs, most of them with a number 
of subsidiary graphs on a larger scale to give increased accuracy. 
Fig. I. Scattered photon energy versus angle, 


"T+ a0(1 — cos6)’ ome 


Each curve gives hy as a function of 0, for a constant initial photon energy hro. 
Fig. II. Recoil energy versus angle, 
2Zaohvo 
1 + 2ao + (1 + ao)? tan? y 


T as a function of ¥, for constant hyp. 
Fig. III. Photon wave length distributior, 


as function of \, for constant Xo. 

Fig. IV. Photon angular distribution, (27)—'eof as function of 0, where 
hv = mc*/) and @ are connected as in Fig. I. Each curve is plotted for constant hyo. 
The last subsidiary graph gives hy as a function of the angle at which oof is 
minimum. 
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Fig. V. Electron angular distribution, 
aof (do, A) (1 + ao)?(1 — cos 6)? i ( 1) 


as a function of ¥, for constant hyo. 
Fig. VI. Photon energy distribution, 


oor” 
F (do, d) 


as a function of hy, for constant hyo. 

Fig. VII. Electron energy distribution, same quantity as in VI as a function of 
T, for constant hyo. 

Fig. VIIIa. Total Compton cross section and effective cross section as functions 
of hyo, Fig. VIIIb. Fraction of incident energy absorbed as function of hyo. 

It is stated that all calculations and curves are accurate to 1 per cent, and that 
the subsidiary graphs are such that interpolated values can be obtained in general 
to 2 per cent accuracy. 

The circular is one of a series of surveys and tabulations of information on 
radiation physics. 

A. E. 


TABLE ERRATA 


242.—R. S. Burtncton, Handbook of Mathematical Tables and Formulas. 3rd 
Edition. Handbook Publishers, Inc., Sandusky, Ohio, 1953. 


I have recently checked, by differentiation, all of the indefinite integrals in 
this edition of the Handbook. The following errors were discovered. They are also 
present in the 2nd edition. 


P. 68, no. 146. In the next to the bottom line of the page, 
for (m + np + n), read (m + np + n)a. 
P. 71, no. 177. In the tan form, insert the restriction: a > 0, c < 0. 
In the tanh— form, insert the restriction: a > 0, c > 0. 


(ad — bc)? (ad — bc)? 
P. 71, no. 178. For + read 
P. 73, no. 195. Insert the restriction: b > 0. 
P. 75, no. 225. Insert the restriction: b > 0. 
P. 75, no. 226. The numerator, —1i, of the coefficient of sin-' U should be 
replaced by: Sgn (d cos ax — csinax), where Sgnz = 1, 
for z >.0, = —1 for z < 0, and =0 for z = 0. 
The expression for U should read: 


v=| 2 + @ + b(ccos ax + dsin ax) |. 
Ve + @|b + ccos ax + dsin ax| 


The final restriction, — * < ax < a, is unnecessary. 
P. 78, no. 258. The restriction should read: b > 0, b > c, cosax > 0. 


42 NOTES 


Most persons using integral tables are aware of errors, which are common to 
many tables of integrals, such as the following. 


The first type is illustrated by the example: f = log x. A correct form for 


this integral would be fe = log |x|. In this latter form, negative values of x 
may be used. 
dx 
The second type is illustrated by the example: f Wns sin'-. This 


form is not valid for a < 0. A form which is valid for all a # 0 may be obtained 
from this form by replacing the a by |a|. 

E. G. H. Comrort 
Ripon College 
Ripon, Wisconsin 


243.—E. JAHNKE & F. Empe, Tables of Functions. Fourth Edition, 1945, New 
York and earlier editions. 
On p. 262, for hi(0.1) = 6.118 read h,(0.1) = 6.342. 
D. A. DE VRIES 
Laboratory of Physics and Meteorology 
Agricultural University 
Duivendaal 2 
Wageningen 
Holland 


Remark: Errors in this volume have been noted earlier in MTAC as follows: 
v. 1, p. 198, 390; v. 2, p. 47, 350; v. 3, p. 41, 314, 364 (review), 423; v. 6, p. 
196 (review 990[L ]), 237. 


240.—See item 2 of the Corrigenda. 


NOTES 
A Conference on Mathematical Tables 


A conference on mathematical tables was held at the Massachusetts Institute 
of Technology on September 15 and 16, 1954, under the leadership of Professor 
P. M. Morse. The following excerpts from his summary will be of interest to 
readers of MTAC. They were written by Professor Morse. 


This conference, under the joint auspices of the National Science Foundation 
and the Massachusetts Institute of Technology, was held to discuss the needs, 
in this country, for Tables of Mathematical Functions, in the light of recent 
developments in high-speed computers. Twenty-eight persons attended the two- 
day sessions and took part in round-table discussions on the general topics : Future 
Need for Tables; What Form Should Tables Have?; What Functions Need 
Tabulating? ; and What Should Be Done About It? 
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General Conclusions 


Still Need for Tables. There was general agreement that the advent of high- 
speed computing equipment changed the task of table making but definitely did 
not remove the need for tables. There are still many scientists and engineers who 
have no more than desk computing machines to help in their calculations, and 
even where electronic computers are easily available many calculations, explora- 
tory and otherwise, still have to be made “by hand.” Here tables of usual type 
and format are still needed. 

Need for a Handbook for the Occasional Computer. There was a general agree- 
ment that an outstanding need is for a ‘‘Computer’s Handbook,” with four- or 
five-place tables of usually encountered functions, together with a discussion of 
their analytic properties and a set of formulas and tables for interpolation and 
other techniques useful to the occasional computer. (This volume was character- 
ized by several as an “‘enlarged and up-to-date JAHNKE-EmpE.”’) The National 
Bureau of Standards for several years has planned such a volume but has not 
been able to afford the five to ten man-years required to prepare, edit, and publish 
it. The present conference strongly recommended that the Bureau produce such 
a volume and suggests that it request the National Science Foundation for finan- 
cial aid to achieve this end. 

Publication of Specialized Tables. It was also agreed that computing machines 
now produce many tables of functions, incidental to the solution of specific prob- 
lems, which would be of use to other workers but which remain in manuscript 
form or in a few copies. The National Research Council publication Mathematical 
Tables and other Aids to Computation (M.T.A.C.) publishes announcements of 
many such tables and keeps copies, which it will loan to interested persons. It 
was agreed by the conference that these tables would be much more useful if they 
were published, even if the publication were by photo-offset or similar process. It 
was recommended that M.T.A.C. extend its services by publishing yearly (or 
biannually) a selection of tables submitted to it. Accuracy of the tables would 
be the responsibility of the producer; the M.T.A.C. editors would exercise the 
sort of supervision over form and content which is done for the usual contribution 
to a scientific journal. The conferees are of the opinion that such a publication 
would eventually pay its way if the printing costs could be held down, but that 
initial cost should be defrayed by the National Science Foundation or by the 
N.R.C. and the Mathematical Societies Tables fund. It is suggested that the 
editorial board of M.T.A.C. (or the National Research Council Mathematical 
Tables Committee) be enlarged to provide assistance in this work; some of the 
additional members should represent tables users as well as table producers. 

“Tables” for Electronic Computers. It was generally agreed that users of digital 
computing machines did not need mathematical tables of the usual form. How- 
ever, those of the conferees familiar with the use of such equipment suggested 
that a collection of algorithms, from which subprograms could be devised for the 
computation of various known functions, when they are needed in a general 
machine program, would be of value. Such formulas or procedures could, perhaps, 
be collected and published by M.T.A.C. or by N.B.S. Computation Laboratory 
from time to time. There also was a present and continuing need for small, high- 
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accuracy tables of key values of functions used by machines for starting series 
expansions or interpolation fermulas which could be incorporated into machine 
programs. 

Additional Suggestions. There seemed to be some demand for a collection of 
tables in punched-card form which can be duplicated or tabulated on request. J 
It was not clear, however, whether the present and proposed collection at Watson 
Laboratories and at U.C.L.A. are likely to be sufficient or not. Further study 
is suggested. 

The need for a continuing Index of Tables was stressed. A new Index will 
shortly be published in England, an up-to-date revision of the Index of Mathe- 
matical Tables, by FLETCHER, MILLER, and ROSENHEAD. It was suggested that § 
M.T.A.C. publish yearly supplements thereafter. The M.T.A.C. editorial board 
will take this suggestion under advisement. 

The rest of the discussion took up the problems of machine production of 
tables, standards of accuracy, format, interpolation techniques, and the kinds 
of functions needing tabulation, reaching conclusions which will be of interest to 
table makers and to organizations supporting table publishing, but which cannot 
now be formulated as specific recommendations. 

Executive Committee. Finally, it was proposed that a committee be formed to 
help implement the recommendations of this report, to cooperate with agencies 
that desire to carry the proposals forward, and to call other conferences if this 
appears desirable. This committee has the following composition : 


Puitire M. Morse, Chairman, Dept. of Physics, Mass. Inst. of Tech. 

M. ABrRAmMowiI7TZ, Div. 11.2, Nat. Bur. of Stand., Washington 25, D. C. 

J. H. Curtiss, 80 Waterman St., Providence 6, R. I. 

R. W. Hammine, Bell Telephone Labs., Murray Hill, N. J. 

D. H. Lenmer, Dept. of Math., Univ. of Cal., Berkeley 4, Cal. 

C. B. Tompkins, Num. Anal. Research, Univ. of Cal., Los Angeles 24, Cal. 
J. W. Tuxey, Fine Hall, Box 703, Princeton, N. J. 


Correspondence concerning this report should be addressed to one of these 
persons. 


International Analogy Computation Meeting 


The Belgian Society of Telecommunications and Electronic Engineers an- 
nounced plans to hold an “International Analogy Computation Meeting” in 
Brussels between September 27 and October 1, 1955. Detailed information may 
be obtained from the Secretary of the Organizing Committee, P. GERMAIN, 
Institut de Physique appliquée, Université libre de Bruxelles, 50, avenue Franklin 
Roosevelt, BRUSSELS (Belgium). 


The Royal Society Depository for Unpublished Mathematical Tables 


The Mathematical Tables Committee of the Royal Society from time to time 
furnishes lists of tables accepted into its Depository of Unpublished Tables, and 
these lists along with short descriptions of the tables are published in the Philo- 
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sophical Magazine and the Journal of the London Mathematical Society. The list 
of tables accepted up to the end of 1953 may be found in the Phil. Mag., s. 7, 
vy. 45, 1937, p. 599-609, and in the London Math. Soc., Jn., v. 29, 1954, p. 504-512. 


Note on Arrangement of Material 


The Editorial Committee is attempting to arrange the contents to make 
reference easier. Generally, the contents will consist of seven sections: 


Papers 

Technical Notes and Short Papers 

Reviews and Descriptions of Tables and Books 
Table Errata 

Notes 

Queries and replies 

Corrigenda. 


In each section material will be grouped according to content approximately 
in order of the Classification of Tables used by MTAC-. 

The most radical change will be the grouping of reviews and notes concerning 
unpublished tables; this grouping is tried as a means of making reference more 
convenient, but it is somewhat influenced by the increasing difficulty of telling 
what has been “‘published” and what is “‘unpublished’’ but widely distributed. 

Another change is the absorption of the portions of the journal devoted to 
automatic computing into the sections listed above. 

No major change in editorial policy is contemplated and none should be de- 
duced from the material in this issue. 


CORRIGENDA 
V. 8, p. 226, 1. —13, for M. W. WiLKEs read M. V. WILKEs. 
for Ch(x, x) = x sin x exp (x — x sin x/sin cosect Add 


read Ch(x, x) = xsin x exp (x — x sin x/sin X) cosect Add. 
. —4, for 70° read 90°. 
. —3, for Ch(x, x) + Ch(x, « — x) = 2 exp (x — xsin x)Ch(x sin x, $7) 
read Ch(x, x) + Ch(x, « — x) = 2 exp (x — xsin x)Ch(x sin x, $x). 
V. 8, p. 227, 1. —15, last column, for .7626582 read .7627582. 
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CLASSIFICATION OF TABLES 
Arithmetical Tables. Mathematical Constants 
Powers 
Logarithms 
Circular Functions 
Hyperbolic and Exponential Functions 
Theory of Numbers 
Higher Algebra 
Numerical Solution of Equations 


Finite Differences. Interpolation 
Summation of Series 

Statistics 

Higher Mathematical Functions 

Integrals 

Interest and Investment 

Actuarial Science 

Engineering 

Astronomy 

Geodesy 

Physics, Geophysics, Crystallography 
Chemistry 

Navigation 

Aerodynamics, Hydrodynamics, Ballistics 
Calculating Machines and Mechanical Computation 


J. 
K. 
L. 
M. 
N. 
P. 
Q. 
R. 
s. 
T. 
U. 
v. 
z. 


CONTENTS 
January 1955 
Coupon Collector’s Test for Random Digits RoBErT E. GREENWOOD 


A Method for the Evaluation of a System of Boolean Algebraic Equations. 


Tables for the Determination of Fundamental Solutions of Equations in the 
Theory of Compressible Fluids 


Analysis of Problem Codes on the MANIAC 
EucGeneE H. Hersst, N. MEtTrRopo.is, & Mark B. WELLS 


Technical Notes and Short Papers 
Wolfram, Vega, and Thiele R. C. ARCHIBALD 


A Note on Approximating Polynomials for Trigonometric Functions. . . 
Continued Fraction Expansion of 2 

JOHN VON NEUMANN & BRYANT TUCKERMAN 


The Values of I'(4) and ['(#) and their Logarithms Accurate to 28 


Modification of a Method for Calculating Inverse Trigonometric Func- 


Reviews and Descriptions of Tables and Books 


GLODEN 1, ONDREJKA 2, PoRTER 3, LINDLEY & MILLER 4, Romie 5, 
STEVENS 6, COHEN, JR., & Woopwarp 7, Iver & Rao 8, ABDEL-Ary 9, 


CADWELL 10, WEILER 11, ROSENBAUM 12, RUBEN 13, TERPSTRA 14, 
Hartitey & Davin 15, BECHHOFER 16, BECHHOFER & SOBEL 17, Dun- 


NETT & SOBEL 18, SENGUPTA 19, SPITz 20, WHITFIELD 21, THOMPSON 22, 
NBS AppLieD MATHEMATICS SERIES 23, RusHTON, RusHTon & LANG 
24, ABRAMOWITZ & RABINOWITZ 25, CLUNIE 26, BERRY 27, BERGHUIS 
28, StarFF oF Comp. Dept. OF MATHEMATISCH CENTRUM 29, NELMs 30. 


A Conference on Mathematical Tables 

International Analogy Computation Meeting 

The Royal Society Depository for Unpublished Mathematical Tables. . 
Note on Arrangement of Material 
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